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MAGNETIC DIPOLE MOMENT DETERMINATION 
BY NEAR-FIELD ANALYSIS 


by 

W. L. Eichhorn 
Goddard Space Flight Center 


CHAPTER I 
INTRODUCTION 


A spacecraft orbiting in the earth’s magnetic field will be subject to forces and torques resulting 
from the interaction between the earth’s field and onboard currents and magnetic matter. Since the 
earth’s field is uniform over the volume of the spacecraft, this interaction may be described solely in 
terms of the dipole moment of the spacecraft. There are two ways to determine this moment. One is 
by direct measurement of the torque acting on the spacecraft in a known magnetic field. The other is 
by a suitable analysis of the spacecraft magnetic field. In cases where torque measurement is impossi- 
ble, one must rely entirely on magnetic measurements. 

This report presents a refined version of the method of dipole moment determination described 
in Reference 1 . It discusses a technique designed to analyze the near-field portion of the spacecraft’s 
magnetic field. The near field is defined as that part of the spacecraft’s field that can not be satisfac- 
torily represented by a centered dipole field. This region begins on the surface of the smallest sphere 
that completely encloses the spacecraft and extends outward until the spacecraft’s field is considered 
“sufficiently dipolar.” This transition from a nondipolar to a dipolar field is not at all well defined 
and frequently occurs at large distances from the spacecraft. 

If, for some reason, measurements cannot be made in the far-field (dipolar) portion of the space- 
craft’s field, data must be taken in the near field. It is here that one lacks a satisfactory data analysis 
technique. The only techniques available for analysis in this portion of the field require data that are 
difficult if not impossible to obtain. Two such techniques are described in References 2 and 3, and 
the difficulties with both of them are obvious. Each requires data to be taken on some closed surface 
surrounding the spacecraft. Whether or not these data can even be obtained depends upon the particu- 
lar spacecraft in question. However, even at best, data collection is tedious and time consuming. Also, 
it would be difficult to devise a data collection system with sufficient flexibility to deal with any type 
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of spacecraft. These are major drawbacks, rendering any near-field analysis technique requiring this 
type of data highly impractical. 

The only points at which the spacecraft’s field can easily be measured lie in a horizontal plane or 
on a vertical line, both of which pass through the center of the spacecraft. If one could calculate the 
spacecraft’s dipole moment from these data, the near-field problem would be solved. Unfortunately, 
an exact determination of the dipole moment from this limited amount of data is, in principle, impos- 
sible. However, it was found that if an assumption were made about the nature of the field, the calcu- 
lations could be done. The spacecraft’s field can be assumed to be adequately approximated by the 
superposition of a certain number of multipole fields. If this assumption is made, a unique solution 
for the dipole moment exists. However, if the assumption is invalid, there will be an error present in 
the calculations. Therefore, this type of near-field analysis must also include a procedure by which 
one can determine the magnitude of this error. If this can be done, one can specify limits within 
which the true dipole moment must lie. Only then is the task complete. 

This report presents a technique which satisfies the above requirements. In the interest of com- 
pleteness, the theory is developed from first principles and carried through to final form. The accuracy 
of the method has been carefully evaluated by computer simulation of many sample cases. The results 
of these tests are discussed along with recommendations for the applications of this technique. 
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CHAPTER I! 


THEORETICAL CONSIDERATIONS 


Magnetically, a spacecraft is a localized source containing macroscopic currents and magnetic 
matter that give rise to a magnetic field. This field obeys Maxwell’s equations everywhere in space. 

The problem is to determine the form of this field in terms of the parameters of the source. There are 
several ways in which one can do this. The formalism which will be used here involves an expansion of 
the field in terms of known vector functions of the angular coordinates, and solution for the coe 1 - 
cients of this expansion by the use of Maxwell’s equations. These vector functions are known as vec- 
tor spherical harmonics” and are most commonly found in the literature of quantum mechanics. They 
are, however, a purely mathematical formalism, useful also for describing classical fields. A discussion 
of them is included in Appendix A and only the directly relevant portions will be employed in the 
main text. It is necessary, however, to start with several preliminary definitions. 

Throughout the text all functions will be referenced to a coordinate system whose origin, 0, is 
located somewhere within the source. The coordinates used will be either spherical polar or Cartesian 
(Figure 1 ). The vector spherical harmonic formalism requires the definition of so-called spherica 
unit vectors.” These vectors will be designated by e +1 , e Q , and e.j, and are defined as 

e +1 = ~ 7^ (u i + iu 2> ’ 

e 0 - U 3 , 

and 

e-l = 7^ (u i"' u 2 ) ’ 

where i = y/ 77 !. These vectors have the following properties: 

(1) The complex conjugate is given by 

e ; = (-iye_. . 

(2) The scalar product satisfies 


(3) Any vector can be expressed in terms of these unit vectors. The spherical components of a 
vector V are given by 

V. = • V , / = ± 1 , 0 . 
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x 3 or z axis 


x 2 or y axis 


or x axis 

r = u 1 cos 0 sin 6 + u 2 sin 0 sin 6 + u 3 cos 6 
6 = u 1 cos 0 cos 6 + u 2 sin 0 cos 0 - u 3 sin 6 
0 = “U 1 cos 0 + u 2 sin 0 

Figure 1— Coordinates of the position vector r in Cartesian and spherical polar coordinate systems. 


If the spherical components of the vector are given, the vector appears as 

+i 

V =]>(-! yvfr 
/=-r 


+i 



/=-! 


In terms of the associated Legendre functions p ™( cos (9), given by the equation 

(- 1 \tn / /i\f +m 

pf'ix) = y- ( i - * 2 r /2 y|) (x 2 - 1 y , ( i > 

ordinary spherical harmonics ^(0,0) can be defined by 
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y . (0, 0 ) = (- 1 ) m \ — y; - ™ lWos ay** . 

/™ v V ^ V 47T 0 + m)! / 


( 2 ) 


The vector spherical harmonics Yj^(0, 0) can now be defined as 

+& +i 

Yf k (6,<t>)= ^ ^ Y k J0,4>)C kl (J,M-m,q) eil , (3) 

m-~k q=~ 1 


where the quantities C kl (J,M\m,q) are known as Clebsch-Gordan coefficients. 

The formulas for the coefficients appropriate to this discussion are given in Appendix A, Table 
A-l . They exhibit the property that C kl - 0 unless k-J ± 1 , / and M = m + q. 

This means that for each value of / except J = 0 there will be three vector spherical harmonics. 
Since the functions form a complete set for scalar functions and the are a complete set of basis 
vectors, the functions in Equation 3 also form a complete set for vector functions in three dimensions. 
Hence, any vector field may be written as a linear combination of the Y^, with coefficients depending 
on the radial coordinate alone. 


Let the steady-state current density distribution within the source be denoted by J(r) and the 
magnetization density by M(r). Let these sources give rise to a magnetic vector potential A(r) and a 
magnetic flux density B(r) related to it by B = V x A. According to the usual conventions, one defines 
the magnetic field strength H(r) by H = B//Uq — M in rationalized MKS units. These functions then sat- 
isfy the following equations: 


and 


J(r), 

inside 

o, 

outside ’ 

M 0 J(r ) + ju Q Vx M(r) , 

inside 

o, 

outside 

ju 0 J(r) + M 0 Vx M(r) , 

inside 

0, 

outside 


V x H(r) = 

V x B(r) = 

V x V x A(r) = 

These equations are valid throughout all space. (Note that J and M are zero outside the source.) 
If A and B are expanded in terms of the vector spherical harmonics, 


°o +/ +i 

A(r) = he d k (j\m,r)Y™. +k (6,<j)) 
/=0 m=~j k=-l 

oo +/ 

= S 1 VM 

/=0 m=~j 
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and 


+/ + 1 


B(r) = S S I 


y-0 m=-j k = - 1 


= I 

j ~ o m--j 


The quantities A^ m and B- m are known as “pure multipole fields”. 

Outside the source, one must require that 

V ‘ B. = V x B. = 0 . 

jm jm 

Therefore (from Appendix A, Equations A-6, A-7, and A-8), 


V * B. 

jm 


and, fromVxB^ = 0, 


}+\ d , / + 2\ 


d' + 


j ( d / - 


2 / + 1 \dr r ) +1 V 2/ + 1 y/f 


2 ) rf:i 


Y. = 0 

jm 


y 2 / + 1 


{[^ G + y 1 ) <i + ^ (I - y 1 ) a] v- 


d j\ . 

— - 1 K y”l, + 

dr r 0 " +1 




'd / + l\ 

7 r + —) d '^,A =°- 


Since the functions and Y ; - m are not zero, their coefficients must be. Hence, 


^ (i + TyK - ’~r) d 'i =0 


so that 


, (± , ill 

\dr r 


) «.* ^ a - y <, - o 


( d j\ , 

— ) d' -0 

<dr r 0 




The solution for Equations 6a and 6b exists only when d' 0 = 0. The solutions for Equations 4 and 
5 exist when 


f d / + 2 n 

<dr r > 
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Therefore, one must require 


d! ~ r~ hl 
u +\ f 


d'-x 


*/-l 


In a region containing no sources, Bj m must therefore have the form 


B /m (r) = — Y™. + ,(0,0) + ¥"..,(0,0) , 


( 7 ) 


r /+2 


where a and b are constants. 

In order to give this form for , the fields must be of the form 

A^- m (r) = d 0 (j,m,r)Y™(6 ,</>) . 


( 8 ) 


This is so because V x Y™. ±1 leads to terms that contain the functions Y which are not found in 
Equation 5. Hence, one is free to take d ±1 (j,m,r ) = 0. 


In order to determine the quantity d Q in Equation 8, the equation 

V*[V*A(r)] = p 0 J(r) + M 0 V x M(r) 
must be solved. From Appendix A, Equations A-3, A-4, and A-5, 


= 7Frr[^S‘7) rf » Y /"-> + ^ (1 


+ 


/ + l\ 

1 ' 7 V m 

0 


and 


Now, 






dr r 


+ (/’+!) 


d i-\\/d , / + 1 


dr 


+ 

r r 


! ) 


d y m 

a Q i. f 


(: 


d j + 2 
Jr r 




d + / + 1 ^ _ _J^_ + _ /(/ + 1 ) 


dr' 


r dr 


as can easily be verified by expansion of the expressions in brackets. Hence, 


V x ( V x A - m ) = ~ 


d 2 2d /(/+1) 
+ 


dr 


2 r dr 


d Q (j,m,r)Y™ = - A(r) J Q (j,m,r)Y”j 


Therefore, 


OO +/ 

V* [ V x A(r) ] = I I - A(r) d Q (],m,r)Y f .(6 ,0) = M 0 J(r) + M 0 V x M(r) . 


7 = 0 m =-/ 


Since 


Ay™ ( 0,0)]*- [Y^;,(0,0)]dS2 = 5 ;/ , 8**. 8, 


(9) 
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the following equation for d Q {j,m,r ) is obtained: 


' d 2 2 d /(/ 

Mr 2 r dr 


d 0 U,m,r) = ~M 0 J Y™ (0; 0')]* [J(V 


) + VxMtr’ddSlIr’) 


- K(j,m,r) . 


(10) 


The function d Q can now be obtained with the use of Green’s functions. Consider the operator A given 
by 


d 1 2 d /(/+ 1) 
A (r) = — + 


dr 


2 r dr 


= ±±( r 2 
r 2 dr \ dr) 


/(/+ 1 ) 


Consider the integral 


/= I g(Af)r 2 dr 


I 

{ g 7r( r2 7/) dr - /(, + 1) f Sfdr • 

‘-'n -'a 


If this is integrated by parts twice, 


I = gr 2 —f 
dr 


t d 

~ fr ~g 
dr 


+ f Jr g ) dr ~ 7 ^ + ^ f ^ dV ' 

0 Jo J 0 


( 11 ) 


If/ and g are required to approach zero faster than 1/r as r approaches infinity, and if they are not 
infinite at r = 0, the first two terms disappear, and 


f 

-'n 


I 


g(Af)r 2 dr = I f(Ag)r 2 dr . 


( 12 ) 


Now, if the solution to the equation 


A(r)g(r,r') 


8(r - r') 


( 13 ) 


were known, the solution to the arbitrary equation 

A(r)f(r)=f'(r) 

could be found as follows: Multiply Equation 14 by r 2 g(r,r') and integrate, so that 


( 14 ) 
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r 2 g(r,r)[A(r)f(r )] dr 


g{r/)f’(r)r 2 dr . 


By Equation 12, this reduces to 


r 2 f{r)[A{r)g{r/)] dr = / f(r)8(r-r)dr 


= /(r') • 


Hence, the solution to Equation 14 is 


•■>=/ 


/O') = / g(r,r')f{r)r 2 dr . 


Therefore, the solution to Equation 1 3 is needed. 


If r = r’. 


£_ + 2d_ _ /(/+ 1 ) 

^2 r dr v 2 \° ] 


g.{r,r) = 0 


The solutions are 


gj(r,r) = 


Ar 1 + A’ r ^ /+1 > forr<r' 
i?V 7 +5r' (/+1) forr>/ 


The boundary conditions on g are, from Equation 11 ,g. approaches zero as r approaches infinity 
and g. approaches infinity as r approaches zero. Therefore, A' = B' = 0 and 


g f (r,r') = 


Ar J for r < r 


Br~ if+1) for r>r 


This can be written as 


r < 

g,(rS) = c — , 


where r > denotes the greater of r and r and r < denotes the lesser of r and r . The constant C can 
now be determined so that the slope of g. has the correct discontinuity at r = r' . 

If one integrates Equation 9 from r = r - e to r = r + e and allows e to approach zero, the value 
for C will have been determined: 
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Lim 

e— >0 


L 


r'+e r r +e 

— r 2 — ^(r,/) dr - + 1 ) / g 7 -(r,r') dr 

- e •'r'-e 


P r'+e 
J t ^ 

r ~e 


5(r - r') dr 


Lim 

e — >0 


? d , 

i d , 


^ ~rg f (r,r ) 

_ jr* g ( r r ) 


dr 1 

dr 1 

r’+e 

t 

r -e J 


= 1 , 
- 1 , 


and 


Lim jL g / rr >) 

e^o df *,<• > •» 

r'+e 


Lim — „ ( r /) 
e ^° dr*i 


, 1\ 2 

-(/+ DC - 

r 


'1\ 2 

=/°7 


Therefore, 


Hence, 


C = - 


2/ + 1 


, 1 r < 
g ^ r) = - 5TT i TTTi 


and the solution to Equation 1 0 can be written as 
d r 


o 


i r r < 

( j,m,r ) = — — / — K{j,m,r)r 2 dr' 

2/ + 1 / -/+! 


> 

OO r j 


' 0 '> 

Outside the source, r > r' always; hence, 


fJL /* 

= — — - / — l[Y%]*-[J(r') +V* M(r')] rfflr' 2 dr' . 

2 ' +I JL ri +1 J " 


1 F 0 f . * 

d 0 (j,m,r) = — — I/'IY %(e'M • (J(r') + V * M(r')] tfF, 

r /+i 2/+1J 77 


where dL is an element of volume defined by r' 2 dr d£l. Now, the solution for A(r) is 


A(r)= S S 1 


/=0 m=-j 


2 / + 1 


77 


7+1 
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where 


D jm = < x 0 y 


>[Y7!(fl',^)']*-[J(r') + V x M(r')] dV . 


n 


The magnetic flux density B(r) and the magnetic field H(r) outside the source are obtainable from the 
relation M 0 H(r) = B(r) = V x A(r). Hence, from Equation A-5, Appendix A, 


+/ 


/=1 m = -j 

which is of the form of Equation 7 with b - 0. 

Using Equation A-10, Appendix A, one may write 
1 


r/ + 2 


ri + 2 


Y™ + i(0,0 ) = [ ( 2/+ 1X/+ l )] * 1/2v 


Lr 


— Y. (0,0) 

y +1 J m 


Hence, 


;=1 m =—j 

= ~ V U(r) . 


/T AAM) 


] jm jm 


+ 1 2/ + 1 r /+l 


The term D . may now be reduced to a more convenient form. Consider the integral 


jm 


f 


[Y “( 6 , 0 )]* Vdfl, 
where V is an arbitrary vector. From Equation A-9, Appendix A, 


where L = -?r x V . Then, 


In the spherical vector notation. 


/> 


/= | [L*Y^ (0,0)1 • V d£2 . 


+ 1 


V* . V = V V*.Vn . . 

1 2 7 2} 


/=- 1 


Hence, 


/ = fia: t r; m >v +l + a; ^ )K. 1 + (Z,; y; )K 0 ] do . 
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Now 


and 


L +i (L x + iL y ) 

e i(t> / a a 

= ~zr [ — - + i cot 0 

R V» 


-iv 

90/ 


1 




e*' 0 / a a 

( — - i cot 6 


/2 \ dd 


L o L z 


1) 

30/ 


so that 


“Z 


30 ’ 


/ 


'+i = I (L +l r. m )v +1 


= - — p + » 17— 
>/2 J LVa® 


icote —J Y J m v^smeded* 


30 


Integrating by parts yields 


- sin 0 f 

2 X 


2 i r 


/ +1 = - — sin 0 j e-t* Y* m (0,0)F +1 </0 


,-2/0 


+ J cos 0 P™(cos 0)F +1 dd 


2 tt 


+ iiy im (M) ^ (sin 0 F+i )£> ~ /0 de d<j> 


J Y fm cos 0 ^ (F+ j e /0 ) ^0 d(j) . 
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Since the first two terms are equal to zero, 


Similarly, 


and 


Therefore, 


Hence, 


/ , = — 1 [y* (d,(p)e~ l<p “ i cot 6 V +1 

+ J2 J [\3<9 90/ 1 

+ ~zj Y jm ( 6 ^ V + 1 ( COS 6 ^ ' COS 6 e ~ 
= -fecMXi;, v +1 )dn. 


/_i dn 

= -J Y Jm <e ^ )(L -l V -0 dn 

! 0 =f V o^l Y ^ e ^ da 

= \ Y ' im W,Wo v ^ dSl - 

1= ~ l \ v -x +L o 

= - jYj m (0»L* • V rffi 

= jY} m (e,<f>)L-vdn. 



Vd£2 = 


1 

\/7U + i ) 



(L* 


sin 0 eft? d(j> 
i(l) )dd d<p 


o)d£l 


V) d£l , 
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so that 


D jm ~ 


M 0 


yj /(/ + 1) 

This can be reduced further; consider 




(0,0)fL • J(r) + L • V x M(r)] dV . 


(r x V)*J=(r,A -x, +(x 0 - x t —V. + (x, — — \J 


2 dx- 3 ax 2 r > ■ r 3 a* A > 3*,p p az " 2 r3 


= Xl (4 /3 - a ! 3 /2 ) +X2 (^ ^ /3 ) + " 3 " 2 - ^ J > 

= r • (V x J) 

= J * (V x r) - V • (r x J) . 

Since V x r = 0, one can write 

L • J = / V • (r x J) . 

Now consider L • (V * M). By the previous equation, 

L • (Vx M) = i V • [rx(VxM)] 

= i V • [V(r • M) ~ M x (V x r ) - ( r . v)M - (M • V)r] 

= / [V 2 (r • M)- V • [(r • V)M] - V • M] . 

Now, 

• a '* 


7)Mi ■ 1 i,(S 



]'k 


. 3 3 3 ' 

5 •/ M . + x , — M 

}k dx. 1 k dx, dx, 1 




Hence, 


Therefore, 


V • M + (r • V )(V • M) . 


L • (V x M) = - / [-V 2 (r • M) + (r • V + 2)V * M] . 




D. =■ 

jm 


y/KT+T) 


rIy Jm * (r x J)] dV- JrJY* m (d,(f))( r • V + 2)(V * M) d/F 


pn 


/ 


+ r'YjmW’Mv (* • m ) c / f , 
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The last integral vanishes upon integration by parts, since V 2 r ] Y*. (0,0) = 0. The second integral can 
also be evaluated by integration by parts. The result is 


/ 

•-'A 


r ] { r * V + 2)( V * M )r 2 dr = -(/+!) I r J V • Mr 2 dr . 


h- 


Hence, 


fX n l 


D jm = 


if 


r'y* m (e^-)V • [r x J + (/' + 1)M] dV . 


\[XT+ D 

The magnetic flux density can now be written as 

•»-* i s wiHwfc 

/= 1 m =-/ 

= -Vl/(r). 

it is now convenient to write this field in terms of real quantities. 


) 


+ M JdV 


fj + 1 


From the definition of Y- m , 


Then, if 


and 


we have 


Y l _ m (e,4>) = (-irY. n (8,v. 


-m Qjm ’ 


+/ / 

^ Q jm = <2 /0 + (2/m + 2^, ) ■ 


m =-/ 


m =+l 


Now, 


^• m (0,0) = 


2 / + 1 U-m)l 
4tt (/ + ra)! 


P/ 77 (cos 0 ) e im<t> = C. P/ 77 e Zi 

/ J*'* J 


im<t> 


Therefore, 


Q. +Q* =C 2 

^ jm jm 


Jpp (cos O’) cos m<f>'f(r') dV PJ 1 (cos 0) cos m0 

+ jpp (cos 0') sin w0 , /(r') dV PJ 1 (cos 0) cos m0 
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One can now write 


oo / 


M 0 /f(cos0) 

B(r) = - ■— V y y (a cos m<p + b. sinra0) 

4tt l— 4 t—* ]m /+ 1 


(15) 


/ = 1 m = 0 


where 


<2. =-(2-o. ) 

m ]0 \i + m 


2 f r ip ” 

)'J 7 


(cos 0) cos m0 V 


and 


2 f r’P V 

OlJ ' 


b- - -2 — — | r'P'" (co'i 6) sin m<t> V 

1 (/ + m , 


" r x J(r) 

L /+1 

r x J(r) 


) 

— 


+ M(r) dV 


/' + 1 


+ M(r) \dV . 


(16a) 


(16b) 


These quantities will be known as multipole moment coefficients. 


Consider now the lowest order term in Equation 15; call this term B x : 

M 0 1 

Bj = ^ [^io^ 0 + (tf n cos 0 + sin (p)P\ ] — 

1 

V [« 10 cos 0 + < 2 n sin 0 cos 0 + £> n sin 0 sin 0] — 


47 r 


M 

47T 


; v hr) ' 


where m is the vector whose components are a lx , b lx , and a 10 . Using Equations 6a and 6b, one can 
now write 


m =|rjvg 


r x J( r ) + M(r) \ ] dV . 


]| 


This vector is defined as the dipole moment of the system. By Gauss’ law, 

J* * [^ r x J(r) + M(r)j| dV- J r £-r x J(r) + M(r) • ds+ J r * J(r) + M(r)J * V J 

surface ' 


volume 


-/[ 


>] 


— r x J(r) + M(r) | dV 

The last form is obtained because the first integral is equal to zero and (A • V)r = A in the second 
integral. Hence, 


m 


=J r x J( r ) + M(r)J dV 
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The significance of this quantity can be understood from the following discussion. Consider an 
element of current density J dV and an element of magnetization M dV in a uniform magnetic induc- 
tion B. The torque on the current element is given by 

dNj = r x [ J( r ) dV * B] 

= r x [J(r)x B] dV , 

and the torque on the magnetic material is given by 

dN M =M(r)xBJF. 

The total torque is given by an integral over the entire source: 


= Jr x [ J(r) x B] x M(r) x B J dV 




J r • J(r )dV+ j I 


J(r)(r • B) dV + B I r • J(r)rfK+ / M(r) x B dV . 


Now, 


Jr • J(r) dV = if 


[ J(r) • V]r 2 dV . 


If this is integrated by parts, 


J'tJCr) • V]r 2 dV = 



surface 


J(r) • ds - 



volume 


[V • J (r) ] 


dV . 


(17) 


The first integral is equal to zero since the current is localized, and the second integral is equal to zero 
since the current is stationary (i.e., V • J = 0). Hence, 




N = / [J(r)(r • B) + M(r) x B] dV 


Consider 


/ 


/= / J(r)(r * B )dV 


This can be rewritten in a more complicated form: 
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1= II J(B • r) - ^ r(B • J) + ^ r(B • J)] dV 


-/[" B ' 


= - / J(B • r) ~ r(B • J)dF + - / r(B • J) + J(B • r)c/F 


/' 


= _ I Jbx ( r x J)dV+ - j r(B • J) + J(B * r) dV 

The second integral can be shown to be equal to zero by the following argument. The zth Cartesian 
component of this integral is 


/ r 


^ [x l (B l J f ) + J l (B / x i )dV= ^ B j f( [x l J,+J t x j ) 
i= i i= i ' 


t x,)dV. 


Now, 

Therefore, 


Vl / = V 


^ Bj f (V/ + •>,•*/) ' (x ‘ 

/-l ^ j=l J 

■ i 4 

/=i ^ 


. V x . + x . V x . ) d. V 

J J * 


J • V x f x . c/ F 


X B i ( J x i x , j - ds ~ j 

j= 1 \ surface volum 


c/s - / X.X. V • J flfF 

* / 

volume 

Both of the above integrals are zero by the same arguments as for Equation 17. Hence, 


/ b " 


N = -— | B x [r x J( r ) ] <iF + J M(r) x BiK 




:/ 1 


- [rx J(r) ] + M(r) J dfjx B 


= m x B . 

Therefore, in a uniform field, the torque acts only on the dipole moment of the system. If the dipole 
moment of the system is known, the torque to be expected in a given field can be computed. 
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CHAPTER III 


ANALYSIS OF THE MAGNETIC FIELD 


Construction of the Near-Field Equations 


It is now necessary to restrict the analysis to the problem at hand. That is, it will be assumed 
from now on that the magnetic field is known only at discrete points in space. The coordinates of the 
zth point are restricted by 

= 0, tt/2, 7T , 

^min ^ r ^ ^*max’ 

and 


0< 0.<2tt . 

Here, r . and r specify the limitations on the radial coordinate and are determined for each speci- 
fic case by physical limitations. The measurement of the magnetic field is usually accomplished by 
rotating the source about the z-axis past a set of fixed magnetometer probes. The axes of these probes 
are aligned with the Cartesian axis of the coordinate system. It therefore makes sense to speak of the 
components of the magnetic field that are seen on the sensors rather than the spherical polar compo- 
nents. Let us now examine these components in more detail and develop the procedure by which the 
quantities of interest (a n , Z? n , and a lQ ) may be determined. 

From Equation 15, the form the components must exhibit can be determined. The procedure is 
to perform first the gradient operation and then, for each value of 6(0, tt/2, i r), make the appropriate 
transformation of unit vectors. These transformations are— 


6=0 

6 =ttI2 

6 = TT 

r = u 3 

r = Uj 

r = ~u. 

6 = Uj 

d = -u 3 

6 = ~u 

0=U 2 

tf> = U 2 

0 = U 2 


Let B.(r, 0, cos 6) designate the component parallel to u ; . for a given probe location. For 6 - 7r/2, 

fJL 00 / 

tfjO-,0,0) y — — (a jm cos m0 + b- m sin m^PJ 1 (0) , (18a) 

l~ 71 r*/^ 2 

/ = 1 m = 0 r 
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n 00 l 

^0 V m 

B 2 (r, 0,0) = “2^ Z, ~~ZZ a jm sin m( & - b jm cos m 4>)P™ (0) , 


/= 1 m =0 


^0 'vp ^ 1 0 

B 3 ^>°) = ^ Z Z '7^ (a / m cos w0 + b f ’‘ sin m4>) Ye p i " (0) ' ( 1 8c) 


j- 1 m = 0 ^ 


For (9-0 and 0 = 7r, the equations are a bit simpler, since at these angles dP m /dd and P m /sin 0 are 
equal to zero unless m = 1 , and the P™ are equal to zero unless m = 0. Therefore, for 0 = 0/ 


JU 00 j 

cos0+6 /l sin ^ 1(1) ’ 


M 00 £* J 

^ fr ’ 0 ' i)= *»*-*,! c ° s «^ 


•p/d) 


47 r r /+2 

7=1 ^ 


(19a) 


(19b) 


NT 7+1 


(19c) 


^ OO 

£ i <'■* > -1 > = cos * + *yi sin *> 1 ) 


(20a) 


. ^°r 1 w 

( n 0 ? 1) - — \ — (a sin 0 - 6 cos 0)— , 

471-4- 7+2 z 1 sin 7T 


/?(-!) 


(20b) 



-(/+!) 
^./+ 2 


pf(-D ■ 


(20c) 


Notice now that a Fourier analysis can create two uncoupled sets of simultaneous equations for 
each value of the index m (except m = 0, for which there is only one set). For a given value of yh, each 
set will involve all the coefficients having that value of m for which / > m. We are interested only in 
the three sets involving the dipole-moment components. The sets involving the coefficients a lx and 
b ll are related only to the fundamental sine and cosine amplitudes of the components B x and B 2 on 
each probe. The set involving a lQ is related to the dc amplitude of the B 3 components. 
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In a Fourier series, the amplitude A m of the sin m<$> term in the function B<4>) is given by 


m 


■if 


27T 


5(0) sin m0 d(j) . 


The amplitude of the cos m0 term is similar, with sin m0 replaced by cos ra0. The dc amplitude is 
given by 


i r 2n 

A ° 2 *J ' 

J o 


5(0) d<j>. 


Therefore, given the measured values of the field components, the Fourier analysis can be per- 
formed easily to yield the systems of equations to be solved. 

Let r i be the radius of the z'th probe location. Let cos 0 z ) be the component of the mag- 

netic field seen on the ;th sensor of the probe at position i. Let ^4 0 (/,/,cos 6.) be the dc amplitude of 
the components 5.(r.,0, cos 0 1 ). Let A x (/,/', cos 0.) be the fundamental (i.e., m - 1) cosine amplitude 
of the component^. (r z -,0, cos 0.). Let A 2 (i,j,cos 6.) be the fundamental sine amplitude of the com- 
ponent 5 ; .(r.,0, cos 0.). These components are given by the following equations: 


i r 2n 

A 0 (i, j , cos 6 f ) = — I 5^,0, cos 0.) d<f> , 

•'n 


o 

277 ' 


A x (i',/,cos 0 


1 f 77 

f) = ~ j 


COS 0y) COS 0 <20 


(21a) 


(21b) 


and 


■ if* 

J o 


AAiJ , cos 0 y ) = - I 5.(r.,0, cos 0 .) sin 0 <20 . 


(21c) 


By use of Equations 18, 19, and 20, these quantities can be related to the various multipole coef- 
ficients. For 0. = 7r/2, the multipole coefficients can be equated as follows: 


oo , 

A i (i ' ] ’ 0)= 'Z i i^ p } (0>a n ■ 


/=. 


A 2 (a,o)=^-^^(0) a/1 , 


(22a) 


(22b) 


/= i r z 


^2 ft 1 > 0) = l Z ^^ (0)6 n > 


(22c) 


r. 

/=1 Z 


21 



^4 ! O', 2,0) = 


(22d) 


and. 


For 0=0, 


and 


For 0 - 7r, 


oo 

- 2 4^’ (0) V . 


/ = 1 r i 


oo 

A o a3,0) = -^~P}(0)a /0 . 


r i 


^ (U ’ 1) = -2 


r i 


oo 

^ 2 ft2,l)= ^ 


1 ^0) 


fl-i 


7 = 1 r i 


/+ 2 sin 0 ;1 




7=1 r i 


/+2 30 


^ i ^/(D 

^ 1 ft2,l) = -> 4-r*,, 

» ~/+2 smO 7 1 


. W 

7=1 r i 




r J + 2 
7=1 r i 


A 


i . _ V 1 3 i 

1 °' 1 ’- I)= 2^^ ( - 1)a / 


j— i r i 


n 


oo 

+ 2 0',2-l)= ]£ 


-P/M) 

— 0. 


^ J./+2 sm 7T 
7=1 r i 


71 


i2(Url)= H^ 1( -‘^ 


7 = 1 


— P^ 72 cos 0 = m cot 0 Pf 72 cos 0 - PJ 1 + 1 cos 0 . 


(22e) 

(23a) 

(23b) 

(23c) 

(23d) 

(23e) 

(24a) 

(24b) 

(24c) 
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(24d) 


and 


A 


oo 

1 (/,2,-l) = - ^ 


P}(- 1) 


v f+2 sin 7 r 
7-1 




^ 0 ft 3 - 1 ) = - y — />“(- 1 )fl /0 , (24e) 

“I” 2 

7 = 1 r i 

It will be shown shortly that the sums in Equation 22 contain only terms for which / is odd, 
whereas the sums of Equations 23 and 24 contain all values of /. However, the symmetry of the func- 
tions Pp can be used to construct equations containing only odd values of /. To do this, the radii of 
the two probes at 0 = 0 and 0 = n must be identical (let it be r A ). The P™ (x) have the following 
properties : 

Pp(-x) = (~iy +m P } m (x) 

and 


Consider now the term 


30 


Pf (-*) = (-l)' +w+1 rr/f (x). 

1 30 7 


1 

2 


1 O', 1 , 1 ) + ^4 j O', 1,-1)] 



-2 

i— 1 


1 9 

r i+ 2 0 ’ 1 3 8 
r A 


-?/(!)[(- 1) + (- i) /+2 i 


Likewise, all terms of the form 


V 1 9 i 

= > rr ^(lH, ■ 

v j+2 30 7 71 

jodd' A 


pA k (i,i,y) + A k {i,i,-V)\ 

can be equated to sums involving only odd values of /. Notice also that since P?(± 1) = 0, 
dPf(+ 0/90 = Pf (+ 1 )/sin 0 and 3 Pf (- 1 )/30 = -Pf (- 1 )/sin n. Hence, 


A j (/, 1 ,k) = -A 2 (i,2,k) 

and 


for k = ± 1 . 
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Therefore, some of Equations 23 and 24 are identical. Let us, then, construct three new equa- 
tions of the same form as those in Equation 22 (i.e.,/ = odd values only) which contain all the in- 
formation present in Equations 23 and 24: 

A X (A)= jU 1 (UJ) + ^ i airl)->l 2 (/,2,l)-^ 2 a2,-l)] , (25a) 

1 

A 2 (A)= -[A 2 (i t l,l) + A 2 (i,l,-1)+A i a2,l)+A 1 (i,2,-1)] , (25b) 

and 

A 0 (A)= j[A 0 (i,3,l) + A 0 (i,3 r l)] ■ (25c) 

These quantities can be equated to the multipole coefficients as follows: 


and 


-A t (A) = 


y > 

Z-r J+2 30 1 
j odd r A 


-A 2 (A) 


I 

j odd 


1 d_ 
j+2 30 

r A 


p}(l)b n , 


(26a) 


(26b) 


+A 0 (A) 


i 

f odd 


y.j+2 


P°(lk 


70 


(26c) 


These equations, along with Equations 22, are the equations that must be solved to determine the di- 
pole moment of the source. They can be placed in a more convenient form if the quantities PJ 1 cos 0 
are replaced with their literal form. For example, it is well known that P?(l ) = 1 for all j. Hence, 
Equation 26c may be rewritten 


A g (A) — 


i 


7+1 


J+2 
/odd r A 


70 


Similar results can be derived for the other values of the associated Legendre functions from the gener- 
ating function for the Legendre polynomials. Since this function is 


F(h,x ) = (1 - 2 hx + h 2 ) 


= 5V*/*) . 

/= o 
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The quantities Pj (0) can now be computed as follows. Consider 

o° 

(1 -x 2 ) 1/2 ~F(h,x)= S h>( 1 -x 2 ) V2 —PAx) 
dx ^ ox 1 

/= o 

h{ 1 -x 2 ) 112 
(1 - 2hx + h 2 ) 3/2 


From the definition of Pj n (x ) 


Equation 27 becomes 


For 0 - 7t/2, x = 0, we have 


/ d\ m 

= (i -x 2 r /2 ( — j p f (x) . 


oo 

^ *'>/<*) 


h( 1 -x 2 ) 1/2 
(1 -2 hx + h 2 ) 312 


oo 

y h>P}(0) = h(l +h 2 f 12 


The right-hand side of this equation can now be expanded in a Taylor series. The result is 


h(l +h 2 )~ 312 =h(\-^h 


2 l.lll 1.1. 111 

+ 2 * 22 ! 2 2 2 3 ! 


Y /!! ( A 0 * 
Z // - l\ \ 2/ 
/0dd \ 2 / 


If this is set equal to the left hand side of Equation 28, 


' U (¥> 


for / odd, 


/’•(0) = o > 


for / even. 


A similar procedure can be followed for dPj (l)/90. The result is 


3 , tfO) l 

-/ , /a) = — 77 =-/(/ + !) 
90 ] sm 0 2 
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These values can now be substituted into Equations 22 and 26 with the index j replaced by the 
index 2k- 1 , where k- 1,2,3,. The resulting equations are: 


^iC/,1,0) - 


t tr 

k= 1 


2k(2k- 1)!! 1 


(/c — 1 ) ! r 2k+l 2k ~ l > 1 


(29a) 


and 


^( 1 , 2 , 0 ) = 


t HT 


(2/c-l)!! 1 


(/c - 1)! r 2^+l 1 » 1 


(29b) 




k(2k - 1 ) 


,2/t+l 


l 2fc-l ,1 ’ 


■4,(U,0) = 


*=1 A 


■ t (-r 


2k(2k - 1 )! ! 1 

(/c- 1)! ’ 


(29c) 


(30a) 


-AA/,2,0) 


oo 

= 2 


/c=l 


1\* _1 (2k- 1)!! 1 


(A— 1)! r 2*+i 2 ^' 1 - 1 ’ 


(30b) 


-AAA) 


A 0 (i,3,0) 


oo 

^2 


k(2k- 1) 


,2k+l 


2/c~l ,1 ’ 


k= 1 \4 


■ 2 a) 

&■= i 


*” 1 (2/c-l)!! 1 


(/c - 1 )! r 2k+l 2k ~ l ’° ’ 


(30c) 


(31a) 


A, 


V 2k 

' (A) ~ L ~Jk7i a2k ~ 1 ’ 0 ' 


k = 1 fa 


(31b) 


Assume now that the magnetic field is measured by N different magnetometer probes. This 
means that the number of equations in each of the above sets is fixed. Since each of these equations 
contains, in theory, an infinite number of terms, no unique solution is possible. One must therefore 
make the assumption that each of these series can be approximated by a finite number of terms. This 
seems reasonable, since the sums are power series in r~ (2k+l K The higher terms can therefore be ex- 
pected to contribute increasingly smaller amounts to the series. 


Specifically, in order to obtain a unique solution, the number of unknowns must be less than or 
equal to the number of equations. Since the lower order terms are expected to be larger than the 
higher order terms, it seems logical to truncate the series rather than to select the various terms to be 
included by some procedure. 
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The number of equations one is able to construct for a given number N p of probes depends upon 
how these probes are deployed. There are basically two configurations that can be considered. One is 
formed by placing all the sensors in the horizontal plane. In this case (known as the + 1 configuration), 
there are 2 N equations for the X- and T-axis moments and N p equations for the Z-axis moments. The 
other configuration (known as the 0 configuration) employs two of the sensors above and below the 
source (i.e., along the Z axis). Here, there are 2 N p - 3 equations for the X- and T-axis moments and 
N - 1 equations for the Z-axis moments. No matter which configuration is used, each equation in a 
set is of the form 


M 

Vt = £ c ifli ’ (32: 

/=1 

where y. is a Fourier coefficient, c t - is a known coefficient, and q- is an unknown quantity propor- 
tional to the zth multipole moment coefficient. 

Solution of the Near-Field Equations 

The problem now is to determine the set of coefficients that best satisfy Equations 32. The cri- 
terion for “best” that will be used here is a least squares one; that is, the multipole coefficients will be 
chosen to minimize the sum of the squares of the deviations of the calculated Fourier coefficients 
from the “measured” ones. 

It is now necessary to define several quantities that will be used in the following discussion. Let 
f(y ) be the probability that if one observes a quantity, that observation will yield the value y for that 
quantity. Usually, the values of y form a continuous set. Then f(y) is defined so that f(y) dy is the 
probability that an observation lies in a range dy centered at y. The expectation value of a function 
q(y) is defined as 


<q>= Jq(y)f(y)dy , 

where the integration extends over all the possible values of y (the integration becomes a summation if 
the allowed values of y are discrete). The population mean, or true value Y of y, is defined as 

Y = <y > 

for a very large number of observations. 

The variance of a quantity is defined as the average of the squares of the deviations from the true 
value. Hence, 

var y = < (y - Y ) 2 > . 

The covariance of two quantities y and z is defined as 


27 



cov (y,z) = <(> - F)(z - Z)> 


■/ 


= I (y~ Y)(z - Z)f(y,z ) dy dz , 


where f{y,z)dy dz is the probability that the observations of y and z lie simultaneously in the range dy 
about y and dz about z. 


In most cases, the probability function /(y) is not known. Also, the number of measurements 
available is usually small. Therefore one must estimate the true values of the measured quantities by 
some method. An unbiased estimate is one whose expectation value is the true value of the quantity. 
Hence, the population mean can be estimated as an average value given by 



i=\ 


Also, the standard deviation of an observation can be estimated from the deviations from the mean v.\ 


Hence, 



y = CW ~ 1 ) var y - {N - 1 )o 2 . 


One therefore estimates the variance o 2 of the observations by the quantity S 2 given by 


The variance of the mean is given by 


N 


^ 0 / - y ) 2 


/=! 


TV- 1 


vary = 


vary 
N ’ 
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Hence, 


S 2 (y) = 


N 

I 


T7\2 


0 / -y) 


i= 1 

N(N-\) 


will provide an estimate of the variance of the mean. 

The above discussion can now be extended to more than one variable. Let y i denote the average 
of a large number of measurements of a variable. Let o f be the standard deviation associated with y { . 

A relationship of the form 

M 

y.-tw (33) 

/-i 

is assumed to exist between the y i and the unknown coefficients q-. The best qj values in a least 
squares sense are needed, along with their associated variances. Assume that there are A known y.. Let 
us call the right-hand side of Equation 33 the theoretical expression for y. and denote its value (for a 
given set of q .) by c r 

The least squares postulate states that the best coefficients minimize the mean of the squares of 
deviations from the means y.. In symbols, 

1=1 


should be minimized. However, the scatter of observations about each of the N points is different 
(of being a measure of this scatter). Therefore, the deviations should be weighted by dividing by of . 
The least squares principle then makes 


i 



a minimum, with w. = o 2 /of and o 2 ~ a constant. 


(34) 


Differentiating Equation 34 with respect to the q. and equating the results to zero yields the M 
so-called normal equations of the system. These can easily be shown to be 


i 

t = 1 


w ^i c ik 




(35) 


These equations can be written in standard matrix notation as 
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(36) 


M 


D, 


= J G kfli > 


/= i 


where the substitutions are obvious from the previous discussion. One method of solution is direct in- 
version of the matrix G. If this inverse is 0, 


M 


1 


0 , * Du = 


M M 

W XT' 


rk k 

k = 1 fc = l / = 1 

M 


Z Z 


1 

7 = 1 


Qr ■ 


Hence, the equations are solved, and the quantities q., given by 


M 


q i = Z */* 5 ’ 

k = 1 i=l 

provide unbiased estimates of the true coefficients 

The variance of these quantities is easily determined, since they are linear functions of the y . : 


(37) 


TV 


° 2 ^ ^ (<ty// 9 T Z ) 2 tf z 2 


z = l 
TV /AT 




M 


i = l \k = 1 




jr W i C ir °i 


\r = 1 


TU Ttf 

■I 

fe = 1 r= 1 


TV 




,7=1 


M M 


X ^ ^ir G kr a ‘ 


k = 1 r = 1 
M 


= 1 */* 6 /*< j2 


fc = l 


= V . 
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The expression for o 2 can now be estimated from the residuals as follows: Consider 


Since 


we can write 

<0 


Now 


and 


Now 


Hence, 


but 


M 


=y r'£ ■ 


/ = ! 


V = 0 


v. - V 

l l 


,') 2 > = 


var v. 


M 


M M 


var vv -2^T w fij cov (y p<?/ ) + ^ co v(q f ,q k ) 

/=1 /=! fc=l 


var v. = oj , 


cov (<?,,<?* ) = <C-? 7 - q, )(<?* - q k )> 


M TV 


M N 


^ X - F r> ^ £ ^p W S C sp O s - r s ) 


\ 


t l=l r=l 
M M N N 


p - 1 5 = 1 


IIIIW, w„, P <0v - };)O s - rp) 


J=1 p = l r=l 5 = 1 


<0-, - *PO s - r,)> = 8 rj Oj . 


M M N 

'ir 


cov 


( v?*> = Z Z Z^ 


c„c , w?cr 2 , 


jr sr sk rv s 


i=l r=l 5=1 


w s =° 2 l°s ■ 
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Therefore, 


Also, 


Hence, 


Now 


M M 


N 


cov 




w s c sr c sk° 


i=l r = 1 
M 


s= 1 


= £ ^/ 6 /jt 0 ' 2 


i = l 


= ^° • 


COV (y,,?,) = <7y,. - r,. )(<?,■ - q , )) 

M TV 

= S S */* c ** w « ^ " y < )(y s " y ^) 


k - 1 s = 1 
M 


-5> 




/c = l 


M M 


M M 


var v 


= <^11 




ij^ik Vjk 

k~l j-1 f—1 k=l 


" 2 S S c v 0 /* c /*° : 


M M 


-•-11 


0., C-.C. 7 O l . 
^]k i] ik 


N 

1 

J=i 


w.r. 


;=1 fc = l 
\ ^ 

: ) = > w i var ^ 

/ “ 

N 

M M 

- I "rt - 

11 

1=1 

7=1 *=1 

M 

= N° 2 -% 

« m ct2 

7 = 1 



N 


W i C ij C ik * ^ 


= (N-M)o 2 . 
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Hence, o 2 can be estimated as the quantity S 2 given by 


S 2 = 


N 

I 

i = 1 


W i V i 


N-M 


Also, the variance of q. is estimated as 


s 2 = 


<p.. JL 

" yr , 2 


N-M 


NT 2 

Z W(V * 


i - 1 


(38) 


A Statistical Test 


Thus far it has tacitly been assumed that any errors in the y. were due to the net effect of a num- 
ber of small disturbances. In this case, the errors will be distributed about the population mean accord- 
ing to the so-called “normal distribution” of statistical theory. This means that the probability of ob- 
taining a value y. will be proportional to 



where is the population mean. 


If x is a linear function of the y. (e.g., a multipole coefficient), X the “true” value of this func- 
tion (corresponding to the F.), and S(x ) the estimated standard deviation of x based on the estimate 


N 

s 2 = ^ 

i= 1 

of o 2 , then it can be shown that the variable 

t = (x - X)/S(x) (39) 

is distributed according to the function 



f{t)dt = 



/ ,2\-(1/2)(v+1) 

v + 7/ 


dt . 


This is often referred to as Student’s distribution, and tables of the variable t can be used to test the 
significance of the departure from an assumed true value. In particular, it can be used to construct con- 
fidence intervals for the value x. For example, for a given number of degrees of freedom v , the prob- 
ability that \t\> t a can be determined from various tables. Let this probability be a. Therefore, it can 
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be said that the probability that the deviation from the population mean is greater than t Q s{x) in mag- 
nitude is a for a given experiment. 

This discussion has provided a method by which the various multipole coefficients and their 
associated uncertainties can be estimated. This method is based upon one basic assumption: The in- 
finite multipole expansion for the Fourier coefficients can be truncated without introducing signifi- 
cant errors. The logical question to ask now is what happens if the fundamental assumption is invalid? 
That is, what happens in the case when one or more of the multipole coefficients, which were assumed 
negligible, actually contributes a significant amount to the Fourier coefficients? The answer is that 
the q. will no longer be good estimates of the actual multipole moment coefficients n.. This is appar- 
ent from the actual expression for the y. : 



M N 

= n i + X X 

k= 1 /=i 

= n, + e, . 

If the basic assumption is untrue, systematic errors can be expected in the q. which cannot be 
eliminated by the above analysis. One can expect the errors e. to change as the order of fit M of the 
system is changed. Specifically, e. can be written as a sum over the multipole coefficients from M + 1 
to °° . Hence, as the value of M increases, more of the coefficients are eliminated, and the error might 
be expected to decrease somewhat. However, this behavior is not definite, and since the error e. can- 
not be separated from the m ., another procedure must be found by which the accuracy of the q. can 
be estimated. 
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CHAPTER IV 


ERROR ANALYSIS 


Introduction 

Consider, for a moment, the currents and magnetic material within the spacecraft. The currents 
are macroscopic currents flowing through coils which permeate the source. Given currents of this type, 
a magnetization density which will produce an identical magnetic field outside the source can, in 
principle, be determined. This magnetization density will not vary appreciably over small distances. 
This is also the case with the magnetization of the magnetic material. The spacecraft source can there- 
fore be divided into macroscopic subvolumes whose magnetizations are fairly uniform. Each of these 
subvolumes can be considered to contain a point dipole source whose magnitude and direction are 
determined by the average magnetization density in that subvolume. By adjusting the size of each sub- 
volume, the magnetic field of this ensemble of point dipoles can be made to coincide with that of the 
actual spacecraft to any desired degree of accuracy. 

Actually, a procedure such as this would be impossible; but the point is clear: The spacecraft can 
be represented by a finite number of point dipole sources, and its field can be written as the sum over 
a finite number of these point dipole fields. It would, therefore, seem that the determination of the 
dipole moment of a point dipole source is of fundamental importance, and that the errors that arise in 
this case are definitely related to those of the general case. 

Let m(k) be the moment of a point dipole located at the point P specified by the position vector 
k. Let the field produced by this dipole be measured by N p probes in either the + 1 or the 0 configura- 
tion, as in Figure 2. The magnetization density M(r) which describes this source is representable by 

M(r) = m(k)S(r - k) . 

This magnetization will give rise to the multipole coefficients of the system as given in Equations 16. 
Upon integration by parts, it will be found that integrals of the form 

J /( r)V • g(r) dV 

volume 

may be written as 

~ J g(r) • V f(r) dV . 

volume 
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ROTATION AXIS 



(a) + 1 configuration. 


ROTATION AXIS 


PROBE 3 



(b) Zero configuration. 


Figure 2— Probe configurations for typical test situations { S p = 2 k/r^ and = r 1 frj) 





( i — 

= (2“5. ft ) m(k) * VrjPf 1 (cos 6) sin ra0 

/U (/ + m)\ 1 

r=k 

(40) 

Similarly, for b jm (k), 

b ■ (k) = 2 — — — m(k) • VrjP™ (cos 6) cos m<t> 

]m (/ + m)! 7 

r= k 

(41) 

The only coefficients of interest in this case are those for which m — 0,1 and / - odd. 
define 

Therefore, 


w ; .(l,k) = a 2/ _ u (k) , 

(42a) 


« ; (2,k) = & 2 /-i,i® ’ 

(42b) 

and 

« ; .(3,k) = fl 2/ . 10 (k) , 

(42c) 


which one can write as 


n.(s, k) = m(k) • h ; .0,k) . 

As before, the actual Fourier coefficients for the 5th axis can be written as 

oo 

^•(5,k)= c..(s)n f (s, k) . 

/= i 

Now, the least squares coefficients q r (s,k ) can be determined as before from Equation 37: 

M(s) N(s) 

q r {s, k) = <f> rp (s) ^ t' ip (s ) w . (5 )y. (5, k ) 

p = i /=i 

°o M(s) N(s) 

■ III 4>rp (s)c ip (s)w i (s)c ij .(s)n f (s, k) 

/= i p - 1 *=i 

oo M(.q at(s> 

= » r (5,k) + 0 rp ( s)w i (s)c ip (s)c if (s)n.(s,k) 

j =M(s )+ 1 p = 1 i = l 
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= n r (sM ) + m(k) • 


« M(s) N(s) 


I II <!>rp ( s ) w i< s ) c ip (s^Oh/s.k) 


f=M(s )+ 1 p~ 1 1 = 1 

= n r (s, k) + m(k) • Q r (s,k) . 

Hence, the error in the coefficient g (s,k) is represented by the term m(k) * Q Cs,k). 


( 43 ) 


Consider now an assembly of point dipoles. The results can be written immediately from a 
summation over all vectors k. Therefore, the least squares coefficient q r (s) is given by 


<I r (s) = ^ q r (s, k) 


where 


= n r (s, k) + m(k) • Q r (s,k) 

k k 

= n r (s) + m(k) • Q r (s,k) 


= n r {s) + e r (s) 


e r (s ) = m(k) • Q r (s,k) , 
k 

if m(k) is discrete, or 

e r (s) =J: m(k) • Q r (s,k) d 3 k , 

if m(k) is a continuous function of the vector k. Now if one could say 

q r (s ) = n r (s ) ± 5 , 

where the magnitude of 5 were known, one would know the limits within which n r (s ) must lie. 
Consider therefore |e r (s)| such that 


Xs)| = 


/' 




m(k) * Q(s, k) d 3 k 


J |m(k) • Q r ($,k)|d 3 k 
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Q r O,k) d 3 k 


< 


/ 


Q?(s, k) / m(k) \d 3 k, 


where Q™(s,k) is the maximum value of the magnitude of Q^fok) for any vector k (within the confines 
of the source). Similarly, for discrete values of m(k) 


k(s)|<Q”fek)^ 

k 

Now, the vector Q r (s,k) and the quantity Q™(s,k ) are independent of the source in question. 
Assume for now that they can be determined by some procedure and hence are known. The error is 
related to the particular source through the terms 

(k) | d 3 k or ^ 
k 

This quantity represents the magnitude of the maximum possible dipole moment of that particular 
source. It is, in general, impossible to calculate. Therefore, some method of estimation must be 
devised. 





The magnitude of this term is dependent upon the magnetic components within the source. If 
they are weak, the term will be small. If they are strong, the term will be large. Such is the case with 
the magnetic field of the object. If the components are weak, the field will be small; if they are strong, 
the field will be large, independent of what the dipole moment is. Also, in the case where 


J^m(k)d 3 k = J |m(k)|d~ 


the dipole moment will usually be large compared with the other moments. Hence, as a fairly good 
approximation, 



m(k) d 3 k 



1/2 


+ B. 


2 pp 




(44) 


where B. pp is the maximum field excursion seen on the / th sensor of the most distant probe during a 
rotation, B . is the maximum value of the field seen on the /th sensor of this probe, and ry is the radius 
of this probe. In the cases where the dipole moment is small, the right-hand side of Equation 44 will 
still be a quantity which will give a measure of the magnitude of the field of the source. The assump- 
tion will therefore be made that 



B . 


(45) 
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in all cases. This assumption is not really critical as long as the estimate of e r (s ) is not so large that it 
is useless, or so small that it is unrealistic. Still, it is an assumption whose validity must be tested fur- 
ther. The discussion of this assumption will be deferred until later. 

Single-Dipole Error Functions 

The form of Q,Xs,k) must now be examined. Evaluate the quantities n } .{s, k); using Equations 40, 
41, and 42, and the Cartesian components of the operator V in the spherical polar coordinate system, 
write the n . as 

= m l [S lf (k) cos 2 (f> k + S' 2 / (k) sin 2 <j) k ] 

+ m 2 t^i/Ck) -5 2 / .(k)] cos sin 4> k + m 3 S 3f (k) cos 4> k , 

Uj( 2,k) = m l [iS 1 ; .(k> -^.(k)] cos <p k sin <f> k 

+ m 2 [S lj .(k) sin 2 (j> k + S 2/ (k) cos 2 <f> k ] + m 3 S 3j .( k) sin <j> k , 

and 

n ji 3,k) = m l S 4/ (k) cos (p k + m 2 S 4f ( k) sin <j> k + m 3 S 5/ (k) , 
where m ; . = the /th component of m(k) and 

k 2j ‘~ 2 T cos e k d 

s !/ (k) = — [ sine * JB 2 /-i( cosfl *:) + Ve p y- i (cos ^> 

^ k 2 i~ 2 p \j-i ( cos 

v ~/(2/-l)L J’ 

k 2 i~ 2 f sin °k a 

^ 3 / -(k) - cos d k ^ 2 /-i( cos 0 *) “ 2j - 1 90 ^ 2 /_i^ cos 

5 4 / (k) = k 2j ~ 2 [(2 j - 1) sin d k P 2j ._ l (cos 0p - cos d k os 0^)3 , 

and 

*^ 5 /(k) = k2i ~ 2 [( 2 / “ 1 ) cos 0 ^ P 2 /-i( cos + sin 0 * ^ 2 /-i^ cos 0^)1 • 

The quantity e r/ .(s) is now defined as 

M(s ) N(s ) 

e r/ S ) = ^ ^ <t> rp (s)w l (s)c tp (s)c l/ (s') . 

P=l 1=1 
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The quantities m(k) • Q r (s, k) can then be written as 

m(k) * Q r (l,k) = m l [E lr (\,\t) cos 2 (j> k + E 2r (\ ,k) sin 2 0^] 

+ m 1 [E lr (l i Y) -/s 2r (l,k)] cos <j) k sin 4> k + m 3 £’ 3r (l,k) cos (f> k ,(46a) 
m(k) • Q r (2,k) = m l [E lr (2,k) - E 2r ( 2,k)3 cos 4> k sin (f> k 

+ m 2 [E lr ( 2,k)sin 2 4> k +E 2r (2, k) cos 2 $ k ] + m 3 E 3r (2,k) sin (j> k , (46b) 

and 

m ( k ) • Q r (3,k) = m l E 4r (3, k) cos <t> k + m 2 E Ar { 3,k) sin <p k + m 3 E 5r { 3,k) . (46c) 

where 


k > = X e n (s > S ,j ( V ■ 

j=M(s )+ 1 

Examination of Equations 29 and 30 reveals that the near-field equations for the x and y axes are 
identical in form. Hence, if the weights w. in Equation 34 are the same for each axis (which they will 
be), e r -{\) - e r j(2 ) for all r and/. Also, since only the errors in the dipole components are of interest, 
all the Ej r {s, k) for which r =£ 1 can be neglected. The dipole errors m(k) • QjCs, k) can now be written 
in terms of only five error functions. These functions are defined as follows: 

E.(k) = E n (\,k) = E fl (2,k), for /= 1,2,3; 

and 

E f (k) = E fl (3,k) , for /=4,5. 

The value of the function Q x (s, k) may be determined for any values of k and 6 k by analyzing 
suitably chosen single-dipole sources. For example, let the source consist of dipole whose moment is 
m (components , m 2 , and m 3 ) located at the point k = ( k , O k ;cp k = 0). Three determinations of the 
least squares dipole-moment components (for three different sets of the ra x , m 2 , and m 3 ) on each axis 
are necessary for a complete definition of the vector Q 3 (s, k) for a particular probe configuration and 
least squares fit. Hence, if m x = 0, m 2 = m 3 = 0, Equations 43 and 46 give 

q x { 1) = m x + m l E l ( k) , 

<7j( 2) = 0, 

and 

q 1 (3) = m 1 E 4 (k) . 

Similarly, if m 2 =£ 0, m x = m 2 = 0, 


( 1 ) = (3) ~ 0 
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and 


q 3 ( 2) = m 2 + m 2 E 2 (k) . 


Also, for m 3 ^ 0, m ] - m 2 - 0, 


q x { 1) = m 3 E 3 ( k) , 


q x { 2 ) = 0 , 

and 

<^(3) = m 3 + m 3 £ 5 (k) . 

Hence, the value of Q^s, k) is known for the vector k. Since the ^-dependence of the function is 
known, Q x ( 5 , k) is really known on a circle of latitude specified by the angle 0 k on a sphere whose 
radius is k. In principle, therefore, the vector Q x (s, k) can be completely determined if this procedure 
is repeated for all values of k. 

Having done this, one must compute 

|Qj(s, k)f =Q 1 fek) • Q x fek) . 


from Equations 46 : 

Qid.k) 


{ [£j(k)] 2 + [^(k)] 2 } cos 2 <p k + [£ 2 (k)] 2 sin 2 <j> k 


Qj(2,k)| 2 = j [^(k)] 2 + t£ 3 (k)] 2 } sin 2 <p k + [£ 2 (k)] 2 cos 2 <!> k , 
Qj(3, k)| 2 = [£ 4 (k)] 2 + [£ 5 (k)] 2 . 


Clearly, 

[<2f(l,k)l 2 < {[^(k)] 2 + [£ 3 (k)] 2 } max + [£ 2 (k>]Lx = 2* ■ 

[Qf(2,k)] 2 < {[^(k)] 2 + [£ 3 (k)] 2 } max + [E 2 {K)} 2 mm =Q l , 

(Gf(3,k)] 2 < {[£ 4 (k)] 2 + [£ 5 (k)] 2 } max =2 3 • 

Therefore, the possible systematic errors {e^s)] in the dipole components can be estimated from 
the expression in Equation 45 for f\m(k.)\d 3 k. Hence, for an estimate of jej(s)| , 


^ 1 (l)| = |c 1 (2)|=C 1 (k)^ 


(47a) 


and 


e x (3) =Q 3 (k)B e . 


(47b) 
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In summary, then, the procedure for dipole-moment determination should proceed as follows. 
Specify a given number N of magnetometer probes deployed in either the zero or the + 1 configura- 
tion about the magnetic source. Then, rotate the source through 360 deg to obtain the three compo- 
nents of the magnetic field at each probe location. Now, Equations 21 will give the appropriate 
Fourier coefficients. If the zero configuration is used, the values of A X (A), A 2 (A), and A Q (A) must be 
calculated from Equations 25. The next step is to solve for the dipole-moment components axis by 
axis. For each axis, first determine the order of fit M(s) desired, and then set up the appropriate near- 
field equations to be solved (Equations 29, 30, and 3.1). 

The least squares analysis now proceeds with the formation for each axis of the matrices D and G, 
as defined in Equations 35 and 36. For the weights w., 1 /of must be used, where a. is the uncertainty 
in the measurement of the values of the magnetic field used to calculate the Fourier coefficient T..* 

The next step is to invert the G matrix and determine the least squares dipole-moment coefficient 
q^s) from Equation 37, and its associated uncertainty S' 2 [<^( 5 )] from Equation 38. Confidence limits 
on the value of q^s) can be determined from the table of ^-values with N(s) - M(s ) degrees of freedom 
[N(s) is the number of equations] and Equation 39. The confidence limits will be given by , 

where a is the probability of a value of t greater than t a . 

Now the possible systematic error in the answer must be determined. First, estimate the maxi- 
mum moment of the source by using Equation 45. Then, determine the value of <2 s (k) appropriate to 
the particular probe configuration used, the sample volume, and the order of fit for the 5th axis. This 
can be done by the rather tedious procedure presented in the last section. Once this is done, the error 
can be estimated from Equations 47a or 47b. 

Obviously, the procedure described would be very tedious if done by hand. However, on a com- 
puter the task is easily accomplished. Let us now turn to the task of evaluating the procedure. 


*These can be taken to be the calibration accuracy of the probe. 
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CHAPTER V 


EVALUATION 


Quantitative T ests 

Thus far, the theory has been developed qualitatively. The equations and procedures have been 
developed without consideration of any experimental limitations that may be present. Also, the fun- 
damental assumption upon which the method is based has been justified by a more or less intuitive 
argument. There is no way to determine how often this assumption can be expected to hold true. 
Therefore, it is necessary that some procedure be devised to test this method in order to determine its 
feasibility as a technique for determination of the dipole moment of a spacecraft. 

The simplest way to do this is to use the near-field method to determine the magnetic moments 
of sources whose actual moments are already known. If the sources selected for analysis are numerous 
and varied, general conclusions can be drawn as to the validity of this technique. Specifically, the test- 
ing program is aimed at answering these questions: Can one expect to determine a reasonably accurate 
value for the dipole moment of an object by this procedure, and is the uncertainty in this number 
small enough to be useful? 

In view of the complexity of the calculations involved in such an analysis, the test program was 
carried out entirely on a computer. That is, a program was written which would generate magnetic 
data for a given set of arbitrarily located point dipoles and then analyze these data by the method pre- 
sented above (see Appendix B). In this way, a large number of different cases could be analyzed in a 
short time. These cases would also be free from any experimental error which could be confused with 
the systematic errors associated with the technique itself. To determine the sensitivity of the technique 
to the experimental errors, provisions were also made in the program to introduce errors into the data 
which would simulate errors encountered in an actual test. 

A modified version of this program was used to determine the error functions ^.(k) on spheres of 
different radii. Hence, the systematic errors could be estimated for each test case by the foregoing 
procedure. 

Test Procedure 

In these tests, it was assumed that the maximum number of magnetometers available would be 
four. The major parameters of interest in a test were the ratio of the source diameter to the radius to 
the closest sensor, and the order of fit of the near-field equations. The parameters varied from test to 
test were the radii of the various probes and the configurations of these probes. 
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The tests were conducted in three sections. The first section was concerned with the determina- 
tion of the moment of a single point dipole located somewhere within a sphere of known diameter. 
These tests were fairly extensive and the results used to compute the error functions E.(k) for a given 
probe configuration, order of fit, and sphere diameter. Also, the order of magnitude of each of the 
lower order multipole coefficients in the near-field expansion could be estimated from these results. 
This is helpful in attempts to estimate the magnitudes of the coefficients of an aggregate of dipoles. 
The results were also helpful in determining the effects of machine error in the computations and the 
behavior of the least squares coefficients as a function of the order of fit. 

The second section of the tests was concerned with multiple magnet sources. Naturally, these 
tests could not be as extensive as the ones above, because of the many possibilities in the construction 
of the sources. Sources containing up to 28 dipoles were analyzed. For a given probe configuration, 

the moments, and their variances, of each of the sources were determined as a function of the order of 
fit. 


In the third section, the same sources as in the second section were analyzed. In this case, how- 
ever, the data were adjusted to simulate real data, and the effect on the calculations was studied. 

In an experimental environment, there will be a minimum observable change in the magnetic field 
because of the limitations of the instrumentation. This threshold value will not be constant but will be 
dependent upon the values of the magnetic field seen on the sensors.* Real data will therefore be sub- 
ject to a minimum uncertainty. Also, the magnetometers have a certain calibration accuracy that must 
be taken into consideration. Errors in the data within these limits were also introduced. 

Obviously, all possible errors could not be simulated. The reason for these tests was to determine 
whether or not errors of the type described previously are detectable from observations of the calcu- 
lated values of the multipole moments and whether or not limits could be placed on these moments 
when this type of error was present. 

Each source was tested using from one to four probes in both the zero and the + 1 configuration. 
This procedure permitted each of the various probe configurations to be evaluated. 


*The magnitude of the fields observed will determine the range setting of the magnetometers. The higher the range the greater the 
threshold value will be. 0 
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CHAPTER VI 


RESULTS 


Accuracy 

The tests show that this method of analysis is indeed effective. They also allow general pre- 
dictions to be made about its capabilities in a given situation. Of primary interest is the accuracy 
of this technique. 

The accuracy of this technique (defined here as the ability to predict the correct values for 
the dipole-moment coefficients) is a function of several variables. These variables can be divided 
into two classes: (I) Variables associated with the magnetic source itself and over which one has no 
control (e.g., strength and location of the magnetic components); and (II) variables associated with 
the technique (e.g., number and location of the probes and the order of fit to the near-field equa- 
tions). The variation in accuracy produced by the Class I variables is the most significant and, 
unfortunately, the most difficult to deal with. That produced by the Class II variables is for the 
most part predictable and is usually much smaller than the variation due to Class I parameters. 
Therefore let us consider Class II first. 

The Class II variables can in turn be divided into two subclasses: (1) physical variables such 
as the number of probes and the configuration of these probes, and (2) the order of fit. There 
appears to be a complicated interrelationship between the accuracy and the physical variables. 

Studies of the error functions E^k) for each probe arrangement possible with from one to four 
probes and for several different sets of probe radii show that there is a definite change in accuracy 
if one of the parameters is changed; however, this change is usually small. The numbers Q l (k ) and 
Q 3 (k) are a measure of the maximum possible error in the measurement of a magnet located on the 
surface of a sphere of radius k (see Figure 2). The characteristics of different probe configurations 
are shown in Table 1; corresponding values of Q 1 and Q 3 are shown in Tables 2 and 3 for different 
size parameters S p (defined as the source diameter divided by the radius to the closest probe) and 
order of fit. 

These results permit one to compare the various probe configurations. In general, a given num- 
ber of probes in the + 1 configuration are more accurate than the same number ol probes in the 
zero configuration (for a certain sample size and order of fit). In addition, the probes provide more 
information* about the magnetic field of the source, and this information permits calculations to be 
carried out to higher order. 


*That is, the number of equations for a given axis is larger than for the zero configuration. 
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Table 1— Selected probe configurations. 


Configuration 

Type 

Number of 
Probes 

Probe Radii (r 1 

=>>/) 


'*11 

r l2 

r l3 

r l4 

1 

+ 1 

4 

1 

2/3 

1/2 

2/5 

2 

+1 

4 

1 

1/2 

1/3 

1/4 

3 

+ 1 

4 

1 

5/6 

5/7 

5/8 

4 

+1 

3 

1 

2/3 

1/2 


5 

+1 

3 

1 

1/2 

1/3 


6 

+ 1 

2 

1 

2/3 



7 

0 

4 

1 

2/3 

1 

1 

8 

0 

4 

1 

1/2 

1 

1 

9 

0 

4 

1 

2/3 

1/2 

1/2 

10 

0 

3 

1 

1 

1 



For probes in the + 1 configuration, it appears that the accuracy decreases somewhat as the 
probe spacing decreases. However, this effect is important only for very close spacing, and in 
general the variation of accuracy is unpredictable. This means that one can assume that the 
probe spacing is not critical. 

The order of magnitude of the change in accuracy as a function of the order of fit is also 
observable in Tables 2 and 3. This change is obviously more drastic than the other changes and 
is definitely predictable. The tables show that as the order of fit is increased, the accuracy im- 
proves. The apparent discrepancy observable in the higher orders of fit is believed due to 
machine errors in the computations. This point will be discussed more fully later. 

The same behavior was observed when multiple dipole sources were tested.* That is, the 
absolute accuracy remained about the same for all probe configurations tested with a particular 
sample and a given order of fit and got increasingly better as the order of fit was increased. 
Again, the discrepancy in accuracy was noted as the order of fit became fairly high, and as be- 
fore, it was assumed to be due to machine error. 

The variation as a function of the size parameter S is also easy to characterize. As the 
size of the sample is increased, the accuracy decreases. This is due to the increase in the magni- 
tude of the multipole moment coefficients as the magnetic material is displaced further from the 
center of the object. This type of variation is really tied in with the Class I variables, which 
are considered later. 

The accuracy as a function of the source parameters is much more difficult to analyze. 

This accuracy is dependent on the values of the various multipole moment coefficients of the 


*The data used in these analyses were considered “perfect”. 


48 



















Table 2— Values of Q l for selected probe configurations. 


Size 

Parameter S p 

Configuration 

Order of Fit 

1 

2 

3 

4 

5 

6 

7 

0.4 

1 

0.13901 

0.0039488 

0.00003915 

0.00001025 

0.00000349 

0.00000995 

0.00000837 


2 

0.12232 

0.0039475 

0.00006992 

0.00000291 

0.00000259 

0.00000574 

0.00000175 


3 

0.12824 

0.0045939 

0.00009143 

0.00001109 

0.00000595 

0.00004668 

0.00056881 


4 

0.13944 

0.0040254 

0.00003694 

0.00001092 

0.00000445 




5 

0.14551 

0.0037651 

0.00007186 

0.00000317 

0.00000402 




6 

0.14089 

0.00419 

0.00001967 






7 

0.081034 

0.018571 

0.00039829 

0.00000672 





8 

0.087116 

0.019120 

0.00036852 

0.00000837 





9 

0.13830 

0.0015225 

0.00011698 

0.00000350 





10 

0.088199 

0.019279 






0.6 

1 

0.29908 

0.01812 

0.00040640 

0.00003516 

0.00000569 

0.0000173 

0.00001552 


2 

0.26348 

0.018058 

0.00071926 

0.00001203 

0.00000423 

0.00000950 

0.00000270 


3 

0.27632 

0.020997 

0.00092398 

0.00003790 

0.00001026 

0.00009265 

0.0011069 


4 

0.29999 

0.018466 

0.00038315 

0.00003621 

0.00000737 




5 

0.31291 

0.017250 

0.00073907 

0.00001258 

0.00000662 




6 

0.3031 

0.0192 

0.00020491 






7 

0.24835 

0.097130 

0.0045043 

0.00003703 





8 

0.26437 

0.099920 

0.0041643 

0.00018160 





9 

0.29747 

0.0067129 

0.0012293 

0.00001845 





10 

0.26738 

0.10074 






0.8 

1 

0.52035 

0.054118 

0.0020232 

0.00030860 

0.00001208 

0.00002730 

0.00002628 


2 

0.45609 

0.054393 

0.0037083 

0.00010345 

0.00000640 

0.00001412 

0.00000369 


3 

0.47795 

0.063327 

0.0048125 

0.00033310 

0.00002863 

0.00017379 

0.0020098 


4 

0.52196 

0.055190 

0.0019150 

0.00031762 

0.00001466 




5 

0.54547 

0.051762 

0.0038109 

0.00010783 

0.00000988 




6 

0.52750 

0.057537 

0.0010723 






7 

0.63746 

0.32796 

0.025616 

0.00039433 





8 

0.67300 

0.33698 

0.023657 

0.0017447 






















Table 2 (continued)— Values of Q 1 for selected probe configurations. 


Size 

Parameter S 

p 

Configuration 

Order of Fit 

1 

2 

3 

4 

5 

6 

7 

0.8 

9 

0.51801 

0.022397 

0.0063067 

0.00016617 





10 

0.68006 

0.33967 






1.0 

1 

0.85543 

0.14429 

0.0079381 

0.0021425 

0.00009264 

0.00004305 

0.0004290 


2 

0.74244 

0.14615 

0.015732 

0.00068123 

0.00001328 

0.00002001 

0.00000468 


3 

0.77575 

0.17022 

0.020620 

0.0023200 

0.00027406 

0.00035247 

0.0036208 


4 

0.85814 

0.14723 

0.0073447 

0.00220801 

0.00009856 




5 

0.89990 

0.13863 

0.016169 

0.00071035 

0.00001591 




6 

0.8676 

0.15383 

0.00373 






7 

1.5286 

0.91619 

0.10295 

0.002561 





8 

1.6045 

0.94099 

0.094992 

0.010326 





9 

0.85274 

0.06508 

0.023189 

0.00088478 





10 

1.6203 

0.94833 






1.2 

1 

1.4622 

0.37326 

0.027578 

0.013383 

0.00080394 

0.00009329 

0.00007073 


2 

1.2464 

0.38442 

0.062195 

0.0039409 

0.00007129 

0.00002966 

0.00000570 


3 

1.2959 

0.44705 

0.082707 

0.014533 

0.0025463 

0.0011433 

0.0062675 


4 

1.4670 

0.38128 

0.024838 

0.013825 

0.00084557 




5 

1.5474 

0.36242 

0.063934 

0.0041143 

0.00005070 




6 

1.4842 

0.40011 

0.010361 






7 

3.7073 

2.4203 

0.34892 

0.012653 





8 

3.8748 

2.4811 

0.32174 

0.046700 





9 

1.4608 

0.18972 

0.079101 

0.0039053 





10 

3.9110 

2.4999' 






1.4 

1 

2.7402 

1.03678 

0.093961 

0.085541 

0.0075137 

0.00099240 

0.00022100 


2 

2.2770 

1.0917 

0.25550 

0.022826 

0.0054193 

0.00006971 

0.00000866 


3 

2.3508 

1.26175 

0.33878 

0.092237 

0.025252 

0.068796 

0.11577 


4 

2.7500 

1.0603 

0.08050 

0.08859 

0.0078824 




5 

2.9211 

1.0220 

0.26269 

0.02386 

0.00033284 




6 

2.7839 

1.1182 

0.031555 
















K> 


l able 3 Values ol Q 3 for selected probe configurations. 


Size 

Configuration 


Order of F 


Size 

Configuration 

Order of Fit 



1 

2 

ra 

Parameter S p 

1 

■ 

2 

3 

0.4 

1 

0.1566 

0.0051880 


1.0 

1 

0.83068 

0.14346 

0.012987 


2 

0.16862 

0.0030402 

0.00002303 


2 

0.87703 

0.081955 

0.0039768 


3 

0.14484 

0.0072891 

0.00025152 


3 

0.74015 

0.21665 

0.046100 


4 

0.1573 

0.0054296 



4 

0.83341 

0.15065 



5 

0.16386 

0.0032509 



5 

0.87833 

0.087407 



6 

0.15876 




6 

0.84315 




7 

0.21510 

0.018625 



7 

2.6890 

0.91858 



8 

0.21749 

0.01904 



8 

2.7140 

0.93646 



9 

0.13605 

0.0034622 



9 

0.73895 

0.11821 



10 

0.21795 




10 

2.7192 



0.6 

1 

0.32077 

0.023574 

0.00084925 

1.2 

1 

1.5463 

0.40008 

0.051806 


2 

0.33434 

0.014003 

0.00024306 


2 

1.6400 

0.21115 

0.012235 


3 

0.29842 

0.032721 

0.0024838 


3 

1.3518 

0.62490 

0.20518 


4 

0.56441 

0.097380 



4 

1.5515 

0.42064 



5 

0.33481 

0.14964 



5 

1.6425 

0.22673 



6 

0.32494 




6 

1.5704 




7 

0.56441 

0.099480 



7 

6.0039 

2.4179 



8 

0.57038 

0.024655 



8 

6.0565 

2.4623 



9 

0.27673 

0.016863 



9 

1.3848 

0.23419 



10 

0.57154 




10 

6.0676 



0.8 

1 

0.50401 

0.064527 

0.0040985 

1.4 

1 

3.0718 

1.1283 

0.20238 


2 

0.52386 

0.039034 

0.0012096 


2 

3.2770 

0.57289 

0.044889 


3 

0.47259 

0.088139 

0.01 1496 


3 

2.6112 

1.8461 

0.89865 


4 

0.50548 

0.67424 



4 

3.0825 

1.1886 



5 

0.52457 

0.041678 



5 

3.2820 

0.61591 



6 

0.51036 




6 

3.1219 




7 

1.2568 

0.32968 



7 

15.396 

6.8312 



8 

1 .2693 

0.33647 



8 

15.525 

6.9503 



9 

0.43033 

0.05086 



9 

2.7694 

0.41841 


. 

10 

1.2718 




10 

15.553 







source under observation; hence, it varies greatly from source to source. As mentioned above, when 
the size of the source is increased, the error usually increases; however, when sources of the same size 
are examined, predictions about the absolute accuracy of the results are in general difficult to make. 

In fact, unless one has more information about the source than its physical size, absolute predictions 
are impossible. That is, it is impossible to say, for a given set of physical parameters, order of fit, and 
source dimensions, that the dipole moment of any source can be calculated to within so many mag- 
netic moment units. This is because the magnitude of the error is dependent upon the strengths of the 
magnetic components of the source: The more “magnetic” the source, the larger the error will be. 

The tests do show, however, a limit (in magnetic moment units) on the magnitude of the error for 
sources of a particular size and magnetization as a function of the order of fit. 

Table 4 shows the maximum errors observed as a function of the order of fit and sample size for 
54 samples containing from 9 to 28 magnetic components (most of which were located at the outer 
edge of the source). These numbers were obtained with the probes in configuration number 1 (see 
Table 1). Since the horizontal and vertical moment calculations are inherently dissimilar, separate 
entries are made for each. The value of the actual dipole moment, M, is shown to the right of the 
observed error 0 . 

This table is not meant to give any general accuracy predictions. However, the magnitudes of the 
magnetic components in the sources are much larger than are usually found in a spacecraft; hence, the 
errors in this table are expected to be larger than those actually seen. If some other source were tested, 
larger errors could very conceivably be obtained. 

Let us now turn to the question of whether or not one can estimate the accuracy of a given cal- 
culation: first, for ideal data (i.e., data for a perfect system) and second, for real data that is subject 
to errors of some type. 

Effect and Estimation of Errors 

In the tests described above, errors in the calculations are either machine errors or systematic 
errors. As one can see from Tables 2 and 3, machine error becomes apparent only at the higher orders 
of fit; this error is due to the rounding off of the numbers calculated by the computer. This effect is 
also illustrated in Table 5, which table presents the variation of the calculated multipole fields* for a 
centered dipole as the order of fit is varied (again for probes in configuration number 1). Note that 
even at its worst, this effect is very small in the dipole moment and can usually be neglected. 

The systematic errors obviously are not negligible; however, they can be estimated by the proce- 
dure described in the text. This was done for each of the sources tested. Table 6 compares the cal- 
culated systematic error limits with the observed errors shown in Table 4. These numbers are represen- 
tative of the results in general in that in every case tested, the observed error was less than the limits 
calculated by the estimate of | e x (^ ) | given by 


*These numbers are the multipole coefficients i) divided by the smallest radius raised to the 2 ] + 1 power and hence are in 

the magnetic field units (nanotesla). 
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Table 4-Largest errors observed in dipole-moment calculations (values in 10" 3 A-m 2 ). 


Size 


Order of Fit 

Parameter 

1 



2 

Z 



% 


5 

6 

7 

S 













p 












— 




Q 

M 

Q 

M 

d 

M 

6 

M 

Q 

M 

d 

M 

0 

M 

0.4 H 

222.3 

2512 

3.6 

144 

0 


0 


0 



■ 

o 


V 

189.0 

972 

3.5 

640 

0 


- 


— 


I ; 

■ 



0.6 H 

477.1 

2512 

18.3 

1604 

0 


0 


0 


0 


0 


V 

444.0 

972 

15.7 

640 

0 


— 


— 






0.8 H 

+795.1 

2512 

58.0 

1604 

2.1 

1850 

0 


0 



■ 

0 


V 

822.0 

972 

41.1 

640 

3.5 

640 









1.0 H 

1146.2 

2512 

148.0 

1604 

6.9 

1850 

2.1 

1604 

0 


0 


0 


V 

1308.0 

972 

111.0 

275 

17.4 

972 

— 


— 


1 




1.2 H 

1499.0 

2512 

353.0 

1604 

18.1 

1500 

11.7 

1604 

0 




0 


V 

1825.0 

972 

287.0 

275 

57.3 

972 

— 


— 






1.4 H 

1855.0 

1500 

881.0 

1604 

57.0 

1172 

67.8 

1604 

5.0 

1604 



0 


V 

2388.0 

275 

761.0 

275 

155.0 

572 

— 


— 


11 





Note: 6 = observed error. 

H = error in horizontal dipole-moment component. 
V ~ error in vertical dipole-moment component. 

M = true value of dipole-moment component. 

0 = errors less than 10“ 3 A-m 2 . 









Table 5 -Calculated x-axis multipole coefficients for a centered dipole for probe configuration 1 (values in nanoteslas). 


Coefficient 

Order of Fit 

1 

2 

3 

4 

5 

6 

7 


80.0 

80.0 

79.999999618 

80.000004821 

80.000039196 

79.993561739 

79.987034380 

a 3 


0.0 

-0.000000507 

0.000019175 

0.000158821 

-0.054924022 

-0.154667296 

a s 



-0.000000202 

0.000024517 

0.000295795 

-0.197767705 

-0.905817656 

a l 




0.000009742 

0.000252243 

-0.338142809 

-2.742426869 

a 9 





0.00079690 

-0.269945145 

-4.308601377 

a \ 1 






-0.081082108 

-3.295361917 

a \ 3 







-0.967441047 




Size 

Parameter 

0.4 

H 


V 

0.6 

H 


V 

0.8 

H 


V 

1.0 

H 


V 

1.2 

H 


V 

1.4 

H 


V 


Table 6— Comparison of the observed errors with the estimated systematic error limits (values in 10 3 A-m 2 ). 


Order of Fit 



e 

E 

e 

E 

e 

E 

6 

E 

6 

E 

e 

E 

6 

E 

222 

522 

3.6 

11.8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

189 

531 

3.5 

22.5 

0 

0 









477 

1137 

18.3 

43.70 

0 

1.8 

0 

0 

0 

0 

0 

0 

0 

0 

444 

1149 

15.7 

103 

0 

3.7 









795 

2005 

58.0 

136 

2.1 

6.6 

0 

1.3 

0 

0 

0 

0 

0 

0 

822 

1918 

41.1 

288 

3.5 

18.3 









1146 

3352 

148 

381 

6.9 

27.2 

2.1 

5.7 

0 

0 

0 

0 

0 

0 

1308 

3359 

111 

266 

17.4 

52.5 









1499 

5821 

353 

1046 

18.1 

131 

11.7 

37.5 

0 

3.8 

0 


0 

0 

1825 

6638 

287 

854 

57.3 

222 









1855 

13730 

881 

3144 

57.0 

378 

67.8 

259 

5.0 

22.8 

0 

2.8 

0 

0 

2388 

7636 

761 

2805 

155 

818 















































|<?,(1)| = |e,(2)| = Q , (k)B e (47a) 

and 

•’ ki(3)| = Q 3 (k)B e . (47b) 

Note that for the lower orders of fit, the estimated errors are quite large. This indicates that in order 
to obtain meaningful results, a given calculation must be made with order 2 or greater. The minimum 
order of fit is dependent on the value of Q 1 or Q 3 for a given sample size and probe configuration. 
Values of Q 1 or Q 3 of the order of 0.01 or less will give acceptable limits. Using these criteria, the 
minimum orders of fit can be determined for the probe configurations given in Tables 2 and 3. 

The results of these tests indicated that the technique is indeed a useful method of dipole-moment 
determination if the data used are almost perfect. In an actual test situation, however, this is not the 
case. Real data will be subject to a much larger roundoff error than computer-generated data. This 
error is due to instrumentation limitations and is dependent on the magnitudes of the magnetic field 
observed at a particular probe location. For testing purposes, it was assumed that the minimum ob- 
servable value of magnetic field is 0.1S nanotesla, where S is a scale factor determined by the maxi- 
mum value of magnetic field H m seen at a particular radius. For a particular value of H m , the scale 
factor S used in the tests can be determined from Table 7. 

The data used in each of the foregoing tests were modified by being rounded off to the nearest 
0. 1*S nanotesla; then they were reanalyzed. As was expected, changes in the resultant calculations were 
observed, but they were in general small when compared with the systematic error already present. 
Table 8 shows the changes observed for the x-axis dipole-moment component for one particular source 
(Sp = 1 -4) measured with probes in configuration number 1 . 

Table 8 shows the important effects of the roundoff process. Notice first of all that the values 
calculated from the truncated data are in general more accurate than the values using perfect data. 

This effect was observed most often in the first and second orders of fit. It is indirectly due to the 


Table 7— Determination of scale factor S. 


Range of H m 
(nanoteslas) 

S 

From 

— 

To 

0 

100 

1 

100 

200 

2 

200 

500 

5 

500 

1000 

10 

1000 

2000 

20 

2000 

5000 

50 

5000 

10,000 

100 
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Table 8— Dipole-moment calculations for a 0.1r l -diameter source 
(actual dipole moment was 1.500 A-m 2 .). 


Order of 
Fit 

Perfect Data (10 3 A-m 2 ) 

Rounded Data (10 3 A-m 2 ) 

Systematic 
Tolerance 
(10“ 3 A-m 2 ) 

Moment 

Statistical 

Tolerance 

Moment 

Statistical 

Tolerance 

1 

3355.1 1 


1897.47 


13730.0 

2 

1221.46 


1467.84 


5194.3 

3 

1555.70 


1536.53 


470.70 

4 

1470.62 

6.30 

1496.34 

1.43 

428.50 

5 

1502.74 

1.06 

1500.80 

1.13 

37.60 

6 

1499.88 

0.03 

1499.48 

0.27 

4.97 

7 

1499.99 

0.01 



1497.97 

2.63 

1.11 


accuracy assigned to the rounded magnetic field values. The uncertainty due to roundoff at each of 
the data points is ±0.055, where S is the scale factor for that particular probe location. This uncer- 
tainty varies from location to location (since the field levels vary). This means that the weights w i of 
each of the near-field equations will be different for each probe location. Since field levels are higher 
on the closer probes, the equations for these probes will be weighted less than those for the farther 
probes. This means that the systematic errors which arise from the closer probes (which are in general 
larger than those from the farther probes) will be correspondingly reduced. In the cases where perfect 
data were used, all the equations were weighted equally, allowing the total systematic error to be 
larger. 

Another effect of the roundoff is a slightly increased error as the order of fit is increased. This is 
noticeable in the table for orders of fit 6 and 7. This effect is believed to be due in part to roundoff 
of the data and in part to machine error. 

Up to this point, the statistical uncertainty in the calculations has been neglected. This was be- 
cause the uncertainty in the perfect data was very small, and the uncertainty in the calculation was 
mostly due to systematic errors. However, when real data are used, the corresponding uncertainty in 
the calculations becomes important. As the order of fit is increased, the systematic uncertainty de- 
creases. In general, this is accompanied by an increase in the statistical uncertainty. In assigning limits 
to the answer, one must choose the larger of the two. This was done for each case tested and it was 
found that the actual error was within the limits predicted unless the order of fit was fairly high (see 
Table 8). 

Note that the uncertainty due to machine error has not been taken into account. The statistical 
uncertainty is due to roundoff of the data, whereas machine errors are due to roundoff of any number 
calculated by the computer. However, as was stated before, the machine errors are important only if 
the order of fit is high. 
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CHAPTER VII 


DETECTION OF ERRORS 


In an actual test situation, the accuracy usually assigned to the data is not the threshold accu 
racy of the instrumentation but the calibration accuracy (which is in general much larger). Errors in 
the data within these limits are to be expected, and their effect on the accuracy of the technique must 
be examined. In order to do this, the data used in the previous tests were modified to simulate linear 
calibration errors and were then reanalyzed. The results show that these errors are indeed significant 
and that their effect increases with the order of fit. However, comparison of the results shows that 
these errors can be.detected from observations of the behavior of the higher order multipole coeffi- 
cients as a function of the order of fit. 

The calibration accuracy assigned to a given data point was, as before, dependent on the maximum 
field level seen at that particular probe location. This accuracy was assumed to be ±S nanotesla, 
where S is the same scale factor as before. The data points were modified by errors of the form 
£)’ = a £) where D is the actual data point, D' is the new data point, and a is the error factor 
(0.99 < a < 1 .01) for a particular probe location and sensor. 

In order to be able to detect the presence of errors, one must know what type of behavior to ex- 
pect of the calculated multipole coefficients as the order of fit is changed. A few tests were run on 
single- and double-dipole sources whose multipole coefficients were easily calculable (from Equations 
40 and 41). Two observations were made from these tests. 

The first was that the actual multipole coefficients of each of the sources decreased significantly 
as the index of the coefficient was increased, even though the size of the source was large (S p - 1 .4). 
In the worst cases, the coefficients decreased by about a factor of 2 as the index was raised by 1. 
Table 9 compares the actual x-axis multipole coefficients* ** (for an x-axis dipole displaced from the 
center of the source by a distance of 0.6^ along the +x axis) with those calculated by near-field 

analysis. 


The second observation was that the changes in the coefficients became smaller and smaller as the 
order of fit increased. These coefficients also steadily approached the true values (until machine error 
became significant). 


*This corresponds to an accuracy of 1 percent of the maximum field level observable on that particular range. 

**Henceforth, whenever “multipole coefficients” are mentioned, it will refer to the actual multipole coefficients (given by Equatio 

16a and b) divided by and expressed in nanotesla. 
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o Table 9 Multipole coefficients for a single displaced dipole, where a. - 2 and A a. is the change in the coefficient as one 

goes from one order of fit to the next. 







These results suggest the behavior one can expect for an arbitrary source. Since the multipole 
coefficients of a source are in fact combinations of the coefficients for single off-centered dipoles, and 
since the latter coefficients decrease rapidly as the index becomes large, the coefficients of a compos- 
ite source can be expected to exhibit the same general behavior. That is, the higher order coefficients 
will in general be much smaller than those at the lower end of the multipole spectrum of the source. 
Also, these higher order coefficients will in general decrease rapidly as the index is increased. It might 
be argued that by careful selection of the locations and magnitudes of the dipoles within the source, 
a particular multipole coefficient can be set to any desired magnitude in relation to the other coeffi- 
cients. This is indeed true; however, the number of parameters that must be specified increases rapidly 
with the index. This means that the probability of a coefficient at the high end of the spectrum (i.e., 
above the fifth coefficient) having a relatively high value will be small if the dipoles are randomly 
placed within the source. The same thing cannot be said of the lower coefficients. Hence, one can 
expect large variations from source to source at the low end of the spectrum. 

The behavior of the calculated coefficients for an arbitrary source is not at all predictable. How- 
ever, in general, large changes can be expected in these coefficients until the order of fit is at least 
equal to the number of the most significant coefficients of the source. After this, the changes in the 
coefficients should be small (until machine error becomes important). If the real coefficients of the 
source decrease as the index increases throughout the multipole spectrum, the changes in the calcu- 
lated coefficients can be expected to decrease as the order of fit is increased. 

The results of many tests on sources whose multipole coefficients were not known in detail seem 
to bear this theory out. Table 10 shows the calculated x-axis coefficients for a source (S p = 1.2) con- 
taining 20 dipoles. Notice the decrease in the magnitude of these coefficients as the index increases. 
Notice also that the magnitudes of the changes in the coefficients in general decrease as the order of 
fit is increased* (up to order of fit 6). These results are typical of each case tested. 

When errors were introduced into the data, either by rounding off the data points or by intro- 
ducing calibration factors, this behavior was changed. The effect of the errors was the same as that 
caused by machine error with the exception that it occurred at a lower order of fit. It was character- 
ized by increasingly large changes in the calculated multipole coefficients as the order of fit was in- 
creased. The magnitudes of the calculated coefficients also became very large. When the data used in 
the previous analyses was rounded off and reanalyzed, the results shown in Table 1 1 were obtained. 

Table 12 gives the results of the analysis of the same data containing calibration errors. Notice 
that the dipole moment values are affected the least by the presence of this error. Table 1 3 compares 
the observed error in the dipole moment with the predicted error limits (either statistical or systematic, 
whichever is larger) for the above cases. Note that in each case, machine error becomes significant in 
the higher orders of fit. Again, these results are typical in that the error limits are in general predicta- 
ble unless the machine error is large. The presence of this error is characterized by large fluctuations 
in the calculated values of the multipole coefficient from one order of fit to the next and by un- 
realistic values for the higher order multipole coefficients (specifically, for Table 12, order of fit 7). 


*The large changes in the dipole and octopole moments from order of fit 2 to 3 are probably due to a fairly large moment (a 5 - 130 
nanotesla). 
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able 10 Multipole coefficients (in nanoteslas) for a 1.2^ -diameter source calculated using perfect data. True dipole coefficient 
1 15.44 nanoteslas, and A a. is the change in the coefficient from one order of fit to the next. 



Coefficient 



A a l 

^3 

Afl 3 

^5 

Aa s 


A a 1 

a g 

A a g 

a n 

Aci l l 

a i 3 


-20.59 67.50 

88.09 


34.62 


-205.41 -235.14 

240.03 29.73 

96.80 -135.60 

38.8 


-230.95 -232.36 


■127.62 -132.69 


-232.68 


■134.51 


-8.25 












<1 

-2.88 10.67 -108.16 

13.55 118.83 

a 3 

Aa 3 

10.40 -207.69 

218.09 

a 5 

Aa s 

-98.24 

a 1 

A a 1 


a 9 

Aa g 


a i l 



A a 


1 1 


1 3 


i 1.2 r x -diameter source calculated using rounded data. True dipole 
^ change in the coefficient from one order of fit to the next. 


of Fit 

4 

5 

6 


7 

-115.57 

-115.63 

-115.40 

- 

■115.87 

7.41 

.06 

.23 

.47 


-233.18 

-233.52 

-231.72 

- 

■237.91 

25.49 

.07 

1.8 

6.19 


-133.62 

-134.52 

-128.10 

- 

164.90 

35.38 

.9 

6.42 

36.8 


-14.89 

-15.89 

-4.77 

- 

116.51 


1.00 

11.12 

11.74 



-0.38 

8.66 

- 

166.91 



9.04 

175.57 




+2.77 

- 

131.45 




134.22 



39.38 
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Table 12-Multipole coefficients for a \ .2r l -diameter source calculated using “uncalibrated” data. True dipole 
coefficient = 1 15.44 nanoteslas, and A a } is the change in the coefficient from one order of fit to the next. 


Coefficient 

Order of Fit 

1 

2 


3 

4 

5 

6 

7 

a l 

-3.03 

10.69 


108.23 

-115.93 

-116.29 

-117.08 

-122.49 

Aa j 


13.72 

118.92 


7.70 

.36 

.79 

5.4 



10.52 

- 

-207.85 

-234.10 

-236.12 

-242.33 

-313.11 

A a 3 



218.37 


26.25 

2.02 

6.21 

70.78 

a s 




-98.37 

-134.81 

-140.01 

-162.21 

-583.13 

Aa $ 





36.44 

5.2 

22.20 

420.92 

a 7 

A a 7 





-15.33 

-21.12 

5.79 

-59.54 

38.42 

-1337.49 

1277.95 

dg 

A a g 






-2.19 

-33.46 

31.27 

-2041.48 

2008.02 

tsT 

<1 

i-H 

rH 

>3 







-9.57 

-1544.59 

1545.02 

a l 3 

-450.36 



Table 13— Observed errors and computed error limits for dipole moments in Tables 10, 11, and 1 2 

(values in nanoteslas). 



Perfect Data (Table 10) 

Rounded Data (Table 1 1 ) 

Uncalibrated Data (Table 12) 

Order 
Of Fit 

Observed 
Error in 

Calculated 

Limits 

Observed 
Error in 

a i 



Calculated 

Limits 

Observed 
Error in 

Calculated 

Limits 

1 

-118.56 

4187 

-112.56 

4187 

-140.52 

4187 

2 

-228.68 

1069 

-126.11 

1069 

-13.36 

1069 

3 

-7.74 

78.96 

-7.28 

78.96 

-8.92 

78.96 

4 

+ 1.06 

38.32 

+0.13 

38.32 

+0.61 

38.32 

5 

-0.23 

2.30 

+0.19 

2.30 

+ 1.07 

2.30 

6 

-0.02 

0.27 

-0.04 

0.27 

+2.06 

1.18* 

7 

-0.01 

0.20 

+0.43 

0.78* 

+8.81 

4.56* 


*StatisticaI error limit (other values indicate estimated systematic errors). 


To recapitulate, the errors associated with this near-field technique vary with the order of 
fit to the near-field equations. They can be divided into three groups: systematic errors, data 

errors, and machine errors. The first of these, the systematic errors, predominate at low orders of fit 
(M ~ 1,2, and 3) and decrease rapidly as the order of fit is increased. The data errors (i.e., those 
caused by imperfect data) are present at every order of fit but tend to become larger as the order of 
fit is increased. 

These two types of errors are similar in that limits can be placed upon them (which were valid in 
every case tested). Machine errors, on the other hand, are caused by roundoff of numbers by the com- 
puter and are very hard to estimate. Fortunately, they are noticeable only at high orders of fit (usually 
M = 6 or 7), and their presence can be detected by careful observation of the results. 
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CHAPTER VIII 


CONCLUSIONS 


This paper has presented a detailed description of a new technique that can be used to determine 
the magnetic dipole moment of a spacecraft or any other magnetic source. The major advantage of 
the technique is that the data required are always easily obtained and easily analyzed (on a computer). 

An extensive testing program was conducted to verify the theory. In all, over 2000 analyses 
were performed, over 500 of which were on composite sources. The results of the tests allow several 
conclusions to be drawn about the technique. 

First, the best results are obtained if all the sensors are located along the x axis (i.e., in the + 1 
configuration). The spacing of the probes is not really critical, but care should be taken to ensure that 
they are neither too close together 0' max A' m i n ^ 2) nor too far apart (^ max /^ m i n ^ 4). 

Second, the magnitude of the errors obtained will decrease as r l (the smallest probe radius) is in- 
creased. However, accurate answers were obtained for sources whose diameters were 1 Ar 1 . If possi- 
ble, r x should be approximately the diameter of the smallest sphere enclosing the source. 

Third, the intermediate orders of fit (M — 3, 4, 5) will, in general, give the answers with the 
smallest uncertainty. For orders of fit lower than this, the systematic uncertainty will, in general, be 
too large to be useful. For orders of fit larger than this, the possibility of large machine errors exists. 
However, higher orders of fit can be used, provided that the changes in the higher coefficients are taken 
into account. Analyses in which the coefficients exhibit large fluctuations or unrealistic values should 
be suspect. 

It should be noted that the z-axis calculations are different from those of the x and y axes in that 
one has much less data to work with. Although the calculations for this axis are just as valid as for the 
x and y axes, they should not be expected to be as accurate. Also, in an actual test situation, care 
should be taken to ensure that there is no zero level uncertainty on the z-axis sensors of the probes be- 
fore data are taken, because this will also introduce large errors. 

In estimating the systematic errors for a given test case, one can use the values of Q and Q pre- 
sented in Tables 2 and 3. To do this, choose the probe configuration and size parameter S^that best 
fit the particular test situation. This will provide only an approximate value for Q l and Q 3 . However, 
since the actual errors were much smaller than the estimated errors for the lower orders of fit, these 
values are not critical. 
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The test cases used to verify this technique were fairly extensive; but it is not assumed that they 
covered all possibilities. However, the sources used were considerably more magnetic than the average 
spacecraft.* Hence, one can assume that the errors obtained were larger (in magnetic moment units) 
than would actually be observed. 

At any rate, the technique was found to work well in every case tested. It is believed that it will 
work well in any case in which one can obtain good data. It is also not unreasonable to assume that 
further study will yield a completely computerized version of this technique. 


*That is, their magnetic components were stronger. 
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Appendix A 


Vector Spherical Harmonics 


Vector spherical harmonics arise when a vector field is transformed by a rotation of its frame of 
reference. This type of transformation is more complicated than that of a scalar field because the com- 
ponents of the vector field are dependent upon the axes of the particular frame of reference. 

The rotation of a frame can be described in terms of an “operator” J = L + S, where L is the an- 
gular momentum operator —ir x V and S is the spin operator, whose components are 

S 2 = ziij x , S 2 = /u 2 x , and S 3 = m 3 x . 


The operators L and S each have their own sets of eigenfunctions. For the operator L, the spherical 
harmonics Y- m satisfy the relations 


and 


L 2 Y. = /(/ + \)Y. 

jm JXJ 7 ]m 


L Y. - mY . 

z jm jm 


For the spin operator S, the vectors e are defined by 


1 

6+1 = ~7f (Ul +/u 2 )) 


and 


= u 3 > 


e -l = ( U 1 -™2> 


such that 


and 


s2e ^ 2e , 

S z% = <? e ? • 


It is convenient to construct eigenfunctions of the operator J similar to 
Since the operators L and S commute with each other, and each satisfies the 


those above for L and S. 
commutation relations 


(L j , L 2 1 iL 3 , 2 1 f 3 ) iL j , (f 3 ? L j ) iL 2 


and 


(S 2 ^ 2 ) ZiS* 3 , etc. , 
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this is easily done by means of the vector addition law. This law states that if Jj and J 2 are two 
angular momentum operators which satisfy the above commutation rules and which commute with 
each other, and if Y im (q) is the angular momentum eigenfunction of the operator J q such that 

and 

eigenfunctions of the operator J = + J 2 can constructed from products of the functions 

y m ( 1)T r (2). This product is an eigenfunction of the operator J z — J lz + J 2 z with eigenvalue 
M = m + r, but is not an eigenfunction of the operator J 2 = Jj + J2 + J| . However, linear combina- 
tions of the products y- m ( 1 )y qr (2) can be used to form simultaneous eigenfunctions of the operator 
J 2 and J z . This eigenfunction will be denoted by G^j q and will be defined by 


G JU, ~ 1 i C jq (.J,M;m,r)y jm (\)y lm V) , (A-l) 

m - — / r-~q 


where the quantities C ia are real numbers known as Clebsch-Gordan coefficients. These coefficients 

/ H 

exhibit the properties that 


C. = 0 unless M = m + r , 

IQ 

and 

J = j + q, j + q ~ 1 > • • • ? |/ ~ Q | • 

The functions G^ jq are eigenfunctions of the operator/ 2 with eigenvalue /(/ + 1) and of the operator 
J z with eigenvalue M. 


Now, eigenfunctions of the operator J = L + S may be constructed using the functions and Y fm . 
(These eigenfunctions are the vector spherical harmonics Y^-j defined as in the text.*) They are 
given by 


+/ +1 

Yf n = ^ '^C n (J,M;m,q)Y jm (8,<P)e q . (A-2) 

m=—j q~— 1 

The values of C.^ can be computed from Table A-l (taken from Reference 4). 

One is now interested in how the vector spherical harmonics behave when operated on by the 
various vector operators (e.g., V * ; V x ). The following equations were obtained from Reference 4, 
and some will be stated without proof. In these equations, the function / is a scalar function of the 
radial coordinate only. 


*The third subscript is equal to 1 because S 2 e q = 2e q = S(S + l)e^; hence, Yjj s becomes Yjty 
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Table A-l— Clebsch-Gordan coefficients C- X (J, M;j, q). 



The above equations can be proven in the following manner. First of all, the expression for the 
quantity V q f(r)Y. m is needed, where V q is the qth spherical component of V. This can be shown to 
be equal to (Reference 4, p. 80) 
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(A-l 1 ) 



Consider now the quantity 

V- C/Y") . 

This is equal to 


+ / +1 +1 

in 


V r Y jm C fl (J,M;m,q) 


m = —j q=— 1 


= I 1 \ Y ,mC n V.M;m,q) 

m=—j q-~ 1 


/ . + | v 1/2 / g *\ +/ +1 


-/ <7 = — 1 



+ 


‘ + 1\ 7 " +1 

T'l I c /1 (/, M; m, q)C jX (j - 1 ,m + q;m, q)Y f _ lm+q 


m ~ — / q = — 1 


However, because m + q =M, the quantities Y. ±} m+q are independent of the summation over m and q. 
Hence, since Clebsch-Gordan coefficients have the property that 


+i +j 


S S = 5 jj’ 8 mm ‘ 


this expression can be reduced to 


V • (/Y § ) - 


/+! \ 1/2 /a 


2/ + 3 


dr r /f Y j+UM 8 J-l,j 


2 i~ 1 


1/2 


/a 7 + 1\ 

W + Vy / w 6 ™./ • 


Substitution of the values / + 1 , / - 1 , and J for / yields Equations A-6, A-7, and A-8, respectively. The 
proofs of the other equations are similar. 
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Appendix B 


Near-Field Computer Program 


A listing of the program used to generate data for the near field of a system of dipoles and also 
used in the near-field analysis of these data is presented in this appendix. 


Let us define B r ,B 0 , and B as the components with respect to an origin 0 of the magnetic field 
B of a point dipole (moment = m) located at the point k (see Figure B-l). The equations for B r , B e , 
and B , are as follows: 

<p 


B r (r,Q,4 > ) 

B e (r,d,<f>) 


3m*(r-k) ( m * r 

0i - 


CO 5 ' 2 

3m • (r ~ k) 


co 


3/2 


CO 


5/2 




m * 0 


co 


3/2 


and 


BJLr,6,4>) 


3m • (r - k) r m 
0 3 


<t> 


CO 


5/2 


CO 


3/2 


where 

\P l = k ( sin d sin Q m cos (0 “ 0 m ) + cos 6 cos Q m ) , 
\l >2 = /c(cos 6 sin Q m cos (0 - 0 m ) - sin 6 cos ) , 
0 3 = £(sin (9 m sin (0 - 0 m ) , 


and 


co = r 2 + k 2 - 2r0j^ . 
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PRO GRAM N RFIELD 

C TO GENERATE NEAR FIELD MAGNETIC DATA 

EQUIVALENCE (ALPHA (1 ), HX (1 ) ) ,(BETA(1),HX(26) ) , < GAMMA ( 1) ,HX(51) ) 

EQUIVALENCE (UA(l),HY(i) ) , (D8 ( 1) , H Y ( 26 ) ) * < DC ( 1 ) , H Y (51) ) 

DIMENSION MX(IOO) , HY(IQO), MZ ( 1 0 0 ? , D ( 1 0 0 ) , TO ( 100) , PD (100 ) 

DIMENSION DVX (4) ,QVY (4) ,DVZ (4) 

DIMENSION R(4), THETA(4), COEF(?3) 

DIMENSION CX(4), SX(4), CY(4), SY(4),DCZ(4> 

DIMENSION CXP ( 4) ,SXD (4 ) , CY0(4), SYD(4), 0CZ0(4) 

DIMENSION qX(8), QY(8),qZ(8), QXJ(8), Q YD ( 8 ) , QZ D ( 8 ) 

DIMENSI ON I TIT LE ( 20 ) .... 

DIMENSION H X ( 3 6 , 4 >, HY(36,4 ), HZ(36,4 ) 

DI MENSION X(36).Y(36), Z(36), PH I ( 36 ) , PH ( 36 ) 

DIMENSION T(10) 

DIMENSION ALPHA ( 25) ,BETA(25) ,GAMHA(25) ,DA(25) ,08(25) ,DC(25) _ 

REAL MUPX ,MUPY ,M UPZ , MUR , MU T , LR , LT , LP ,MUP 

REAL MX,HY , HZ , ML Y , MLZ , ML X , ML ♦ MU RX , MURY , M UR Z , MUT X , MUT Y , MUT I 

REAL MAX, MIN 


C ZERO OUT ARRAYS 

DO 450 I = 1,4 

CX(I) = C Y ( I ) = SX ( I) = S Y ( I ) = D C Z ( I ) = 0, 

CXO(I) = SXD(I) = CYO(I) = SYO(I) = DCZD(I) = 1. 

ALPHA(I ) = BET A (I) = GAMMA(I) = 0, 

DVX ( I ) = DVY(I) = OV Z( I) = 1, 

DA (I) = ) B ( I ) = DC (I) = 0. 

IT = I + 4 

ax (I) - Q Y ( I ) = QZ ( I ) = 0. 

ax (II) = QY (II) = QZ(II) = 0 . 

QXD(II) = QYD(II) = Q7D (II) = 1. 

QXD(I) = QYD(I) = QZD ( I ) =_ _.i . 

H X ( I ) = H Y ( I ) = HZ<T> = 0, 

450 CONTINUE 

T ( i) - 6. 314 

T (3) =2, 353 ___„ 

T(2>=2.9?0 

T (4) = 2.132 

T (5) =2. 01 5 

T (6) =1.94 3 

T ( Z > =1 . 8° 5 

_ T ( 8 ) = 1 . 86 0 

T(9>=1.833 

T ( 10) = 1.812 

' 3 FORMAT (6 r i0 . 0) 

: 1 FORMAT ( 8E 10.0) 

, 2 FORMAT ( 10 15) 

till READ (60,2) NUM,ICON,NORQER 

I F ( NUM ) 4 3 3,430,350 

< 350 READ (60,1) ( R(I),I = t,NUM> 

* 428 READ (60,?) NO 

IF ( NO) 11.11,1111,429 

4 2 9 RE A D ( 6 0 , 3 ) ( MX ( I ) , H Y ( I ) ,M Z ( I ) ,0(1) , TO ( I ) , PD ( I ) , I = 1 9 NO ) 

JO = 0 _ 

DO 4061 I = 1 , NO 

1 Li R ( 1)-D( I ) ) 40 6 2_, 4 0 b 2,40 61 

4 0 62 " W'p I TE (fei, 40 6 3 ) I 
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40 8 3 FORMAT { ?QX, 1 2H ERROR IN 0 ( > 1 3 , 2 H > ,/) 

' 

40 oi CONTINUE 

IF(JQ) 4^65,4065,423 
4065 CONTINUE 

DO 30 0 r = 1,1 UH 

30 n THfcTA(I) - 90. 

IF (ICON) 301,301,302 
311 THETA (NUN) - 180. 

THtTA (MU’ • - 1) - U. 

30? CONTINUE 

5 HO 10 1 = 1. ,NUM 

10 THETA (I) =THETA (I) *.017 45 32 7 
HO 625 T ' 1, NO 
T D ( I ) = TO Cl) *.0 1 7453? 7 

PD(I) = P0(I) *.01745327 

625 CONTINUE 

NTERMS = NTERM2 = 1 
NP = NUM 
LQ = 0 

JJ = 0 „ 

DO 521 L = 1,36 

521 p H ( L ) 5 ~ L O A T ( L - 1 ) * 1 0 . 

no 00 J=t,NUH 
PHI =-.174 5 32 7 
DO 601 L = 1,36 

601 X ( L) = Y ( L ) = _2 ( L_>_ = _n_._ 

LQ = 1 

DO 100 K= 1 , 56 

PHI=PHIf. 1745327 

PH I K= PH.I 

H P = H T = H P -- 0 . 

RJ = P ( J ) 

S I M T J = S I N (THETA (J) ) 

COSTJ = COS ( THET A ( J) ) 

DO 8 0 1 = 1 , NO 

C CALCULATE MAGNETIC FIELD A T R J , T H F T A ( J ) ,PHI(K) DUE TO PGLfc 

c S = S I NT J * sTnTt ! D ( I) ) *C O ST PH I K - P O (I) ) 

ALPH =CS' COSTJ*COS ( TD( I) ) 

BET =0050 J* SIN (T 0(1) ) *COS ( PHIK-PO ( I) ) -S INT J *COS ( TO ( I ) ) 

G A H M =SIN (TO ( I) ) »SIN ( p HIK~PD ( I) ) 

Q = RJ* RJ 4 - 0 ( I ) *0 ( I ) -2 . D ( I ) * ALPH 

CS=RJ»S INTJ *COS(PHIK ).- SI N ( TO ( I ) ) »C OS ( PO < I) ) »D ( I) 

ML X = MX (I) *CS 

MLY = MY (I) * (RJ*SINT J* SIN (PH IK) -D ( I) *SIN ( TD ( I > ) *SIN (PO ( I ) ) > 
ML Z = M 7 ( T) * ( R J “ COS T J - 6 ( I ) * C 0 S ( TD (I) ) ) 

ML=MLX»ML Y + NLZ 

CS=3.*ML/ (SOR1 (Gl>*Q*Q> 

_ _MURX = MX (I ) *SINTJ*CO'S (PHIK) 

MU RY = MY ( T ) *SINTJ*SIN (PHIO 
MUR7 = M2 (I > *COSTJ 

~ mutx=mx a T» costj*c 6 s (PHIK ) 

MUTY=MY (I ) »COSTJ»SIN (PHIK ) 

MUTZ=-MZ f I) *SINT J 

MUP_X = -MX(I) ^ SIN (PHIK) _ 

M UP Y = M Y (I ) * C O S ( P H I K > 
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MUPZ=Q. 

MUR=MUPX^ MURY+MURZ 

MUT=MUTXf MUTY4-MUTZ 

MU P =MUPX+ MU PY+ MU PZ 

Ql=l . / (SORT (Q) *Q> 

lr=mur*q_{ 

_LTfMUT*Q1. 

LP-MUP^Q1 

C MAGNETIC FIELDS 

BR = CSMRJ-D (I) *ALPH >-LR 
BT= r CS v O( I) ^BET -LT 
BP = OS *D ( I ) * G A MM - L P 

C TO T_AL FI- L 0 

HR=HR+RR 

ai.= HT+.BT.... . 

HP=HP+BP 

30 CONTINUE 

0 1 FF=T HE! A ( J) -1.5/03 

IF(DIFF) 20,30,40 

PC IF <A8S(DIFF) - . 01) 30,3 0,25 

25 Z ( K) = - H Rf_l Q 0 00 0 . 

X < K ) =H ! *100000 . 

Y_(KJL _* "H P_* 1 0 0 0 0 0_* 

GO TO 36! 

30 X ( K ) =H^*10Q0 00 . 

Z ( K) =H T *1 0 0 U 0 0 . 

Y ( K) s-HP* 100 000. 

GO TO 260 

4 0 IF ( A 8S (_D T FF ) - . 0 1 ) 3 0 ,_3 0 , 4 5 

45 Z(K) =HR* 10 00 00. 

X ( K) =-HT»1000Q0. 

Y (K) =-HP*10000Q. 

360 CONTIN UE 

101 FORMAT ( 10X , 3 ( F 2 0 . 1 0 , 2X ) , / > 

_ 1 U 2 FORMAT ( IPX, 1 5_) 

H X ( K » J ) - X ( K) 

H Y ( K , J ) = Y ( K) 

HZ(K,J) - Z(K) 

100 CONTINUE 

C HAVE MAGNETIC FIELD 

C NOW CALL FOU RIER AN ALYS IS SUBROUTINES 


372 FORMAT (/» 20X,L9HCQEF FOR PROBE NO. ,11) 
C 


c 

CALL FOURIER (NOROER, 

X ,COEF) 



CX(J) = C OEF ( 2 ) 
SX(J) = C OEF ( 3 ) 

CALL FOURIER (NORDER, 
S Y { J ) = C0EFC3) 

Y , 0 O £ F ) 



C Y ( J ) = C OEF (2 ) 




CALL FOURIER (NORDER , 

Z , COE F) 



O C 7 ( J ) = COEF(i) 




330 CONTINUE 





IF (N P- J ) 5 0 4,504 ,50 3 
503 IF ( 4- J) 3 0 4,5 0 4,90 


504 JJ = JJ+t 
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O O D O 0|0 


HP - NP-4 
LO = 0 

ISTART ■= (JJ-i) f J «- 1 
I S T CP = J J * 4 " 

WRITE (61 , 155) (I, I -I ST ART, I STOP) 

155 FORMAT ( 1. HI , 4 < 16 X , 5HPR0BE , I? , 1 D X) , // ,6H 'THETA, 6X, 1HX,9X , 1HY ,9X , ihi 
1 , 3 < 12X , l^JX , 9 X , 1H Y , 9X , 1HZ >_/ ) 

no 2245 K ' = 1,36 

WRTTF (61,156) PH (K) , (HX (< , L ) , HY (K , L ) ,H Z (K , L> , L = 1,4) 

2245 CONTINUE 

156 FORMAT (lY,F5.L,3(lX,F9.3),3(4X,f L 9.3,lX,F9.3,lX,F9.3)) 

IF (ICON) 4 0 20,4020,4030 

4020 NJ = NUM - 2 
GO TO 4 031 
4030 NJ = NUM 

40 31 C A L L MAX MI N ( MA X , M IN , 36 , HX ( 1 , N J ) ) 

R3 = R(NJ)*R(NJ)*R(NJ)*.a0001 

HX = .25* (MAX-MIN) *R3 

CALL MAXMIN (MAX,MIN,36,HY (i,NJ) ) 

BY = .5*< MAX-MIN) *R3 

CALL MAXM IN (MAX, MIN, 3b, HZ ( 1,NJ) ) 

' 81 = APS (MAX) 

92 = ABS(MIN) 

IF (91- B 2) 40/0,4070,4080 

40 5 0 97 - MAX* 0 3 
GO TO 40/1 

4070 97 - MIN* R 3 

4071 CONTINUE 

9= SQ°! OX * B X + BY* 9 Y+BZ *87) 

“ 5 01 OO 50 2 I = 1 ,4 

OO 50 2 L = 1,3.6. _ 

H X ( L , I ) - HY ( L , I ) = HZ (L, I) - 0. 

5j? CONTINUE 

2250 CONTINUE 

C 

90 CONTINUE 

AT THIS p O I N T ALL CALCULATIONS ARE COMPLETE 
FPqm HFR- ON IS MAINLY OUTPUT- ^LOTTING AND PRINTING 


PRINT SECTION 


“ CHANG E PH I ~ FRO M R ADI A NS TO O C G R . t L S 

pHi(i)_=a. _ 

OO 150 1-2,36 

19 0 PHI (I)=P p I(I-l)+li. 

IF (17. EQ, 1) 1112,1113 

1112 WRIT E (61, 151) 

151“ FORM AT ( ////////// ,2U X , 45H MAGNETIC FIELD C A LCUL A T I O NS- -PROBE LOCATx 
_ 1 O N S_ ) 

DO 161 1 ,NUM 

THETA ( I ) - THETACI)/ . 01745327 

161 CONTINUE 

_ WRITE (61, 152) ( I , R ( I ) , THETA (I) , 1=1, NUM) 

1 5 2 FORM AT ( 40 X , 6HPR0QE ( , 1 1 , 2 0 H ) IS LOCATED AT R = , FI 1 ,4,136 CM . THLTa 
1 = ,F.11.4,lX_,7HDEGPt:ES /) 

1113 WRITE (61, 153) NO 
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153 FORMAT ( // , 2 OX, 18HFIELO PRODUCED J3Y ,12,2 2H DI°OLfcS DETERMINED BY ) 
DC 162 T = 1, NO 

TO ( I ) = T0(I)/«Q1745327 

PO(I) = PD (I) /. 01745327 

162 CONTINU E 

WRITE (61 , 154) { I ,MX ( I ) , I ,MY ( I ) , I ,H7 ( I ) , I , 0 ( I ) ,1 , TO ( ii , I , PO ( I > , 1 = 1 , 

1ND) 

154 FORMAT ( lx , 3HMX(,I2,2H)=,Fll«4,lX,3HMY<,I2,2H>=,Fli.4,iX,3HMZ(,I2, 
12H)=,F11.4,1X,2HD(,I2,2H)=F11.4,1X,3HT D(,X2,2H )=,F11.4,1X,3 HP0< , 12 

2,2H)=,Fli .4) 

DN= 0 ♦ _ _ 

0 E - 0 • 

DP=0 . 

DO 1010 1=1, NO 

DN= DN + MX ( I ) 

D£ -DE-MY ( I) 

1010 OD=DD-MZ ( I ) 

W R, I T E < 6 1 , 10 11) DN,DE,nD 

WPITECBl, 4032) B _ 

4032 FORMAT<20 X, 2 7HMAGNE T IZ AT ION ESTIMATION - ,F?0.iG,//) 

1011 FORMAT <9 ( /> ,45X, 20 HTOTAL DIPO LE MOMENTS , / / , 5 0 X , 7H NOR TH = , 

1F11.4, / ,5QX, 6H2AST = ,F1I*4,/, FOX, 7HDOWN = , F11.4,////) 

C PLACE F.C«*S IN PROPER FO RM FOR MATR IX 

IF (ICON) 401,401,402 

401 N = NUM-? 

C X A 3 = .25* (CX (NUM-1) +CX(NUM) +SY (NUM-1) + S Y (NJM) ) 

S X A 3 = . 2 5* (SX ( NU M-1 ) + SX CNUM ) - C Y( NUM-1 ) -C YJNUM ) ) 

OCZA3 = . 5* (OCZ (NUM-1) + OCZ (NUM)) 

CX A 30 - ,5 

SXA3D = * 5 

0CZA3D = . .7U7 - - 

GO TO 403 

_4 0 2 N - NUM .... 

403 DC 40 4 I = 1 , N 

NN ~ I + N . ... 

QX ( I ) = G X ( I ) 

GX( NN) = - S Y ( I ) 

Q Y ( I ) = SX(I) 

QY(N N) = CY (I) . 

Q Z ( I > = 0 C Z ( I ) 

4 34 CONTINUE 

IF (ICON) 405,405,406 

405 NM = 2* (MUM-2) 

Q X ( NM+ 1 ) = -CX A3 

Q XO (NM+1) =C X A3 Q - 

QY (NM + 1 ) - -SXA3 

QYO (NM» 1) - S X A 3 D 

Q7 (N+l) -- DC7A3 

070 (N + l ) = ... DC Z A 3D 

NFQX - NM+1 

NEQ.7 = NUM-1 

GO TO 807 

40 6 NFQX = 2* NU M 

NEQ7 = NUM 

C H OW EVERYTHI NG IS IN THE COR RECT FO RM FOR SOLN RY MA TRI XX 

807 CONTINUE 
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WRITE <6_ i, _ _ 

4)7 FORMAT (20 X,2HR , 13X, 3HQX , 12X, 3HQXD , 12X , 3HQY , 12X, 

1 3HQYQ , 12X , 3HQZ , 1 2X, 3HQZU , //) 

911 FORMAT ( 15X, 7 (Fly. 7 , IX j , / ) 

9 1 2 FORMAT ( 1 5X_5__5 (F15. 7 , IX),/) 

W PIT F (61, 91 1 ) ( R (i ) , Q X ( I ) , 1X0(1) , OY ( I ) ,GYD (I ) ,Q7. ( I) ,QZO ( I ) , 

i I =1,N> _ 

DC 913 I = 1,N 

I ? = I 4- N 

WRITE (61, 912) R(I) , QX < I?) , QXD (I2),QY(I?) ,QYD(I2) 

9 13 CONT INUE 

IF (ICON) 41 3, 41 0,411 

41 n KL = _2*NUM -3 

WRITE (61* 911) (RCN + 1) ,QX(KL > ,QXO (KL) , GY ( KL ) , QYO ( KL ) ,GZ <N + i) , 

1 Q7D (N+ 1) ) 

411 CONTINUE “ “ ~ ’ 

4 26 CALL M A T V I X ( IC ON , 1 , NUH ,NEQX , NT FFMS , R , QX , QXD , ALPHA , DA) 

42 0 CALL MATR IX ( IC ON , ! ,~NU- , NFG X , NT E P IS * R » Q Y , QYO, BET A , DB) 


3 2 3 FORMAT ( 

SOX, 

6HXAXIS 

y 

/ ) 

804 FORMAT ( 

SOX, 

bHYAXIS 

Vn 


ft 0 5 FORMAT ( 

POX, 

GHZ AXIS 

9 

//) 


CALL M A TR I X < ICON , Q , NUM , NEQ7 , NTERM? , R , QZ , QZ D , GAMM A ,^(T) 

WRITE (61, 27 ) 

DO 50 'I = i , N T ER M S ' " " 

J = 2 * I ~ 1 

WRITE (61, 21) (J , ALPHA ( I ) , D A (I) , J , IF T A ( I ) ,09(1) ,J, GAMMA (I) , 'DC (I) > 

21 FORMAT ( 4 H A(,I2,6H,1) = , F 1ft . 9 , 1 X , 2H 4 - , F 1 0 . 5 , 

1 4 H B(,I2,6H,1) = ,F18.9,1X,2H+-,F10.5, 

1 4 H A ( , 1 2 ,6H , 0 ) - ,Fift. 9 ,1 X , 2H + - ,F10 . 5) 

>0 CONTINUE 

27 FORMAT (30X, 32HC ALCULA TEO MULT I ROLF COEFICIENTS , //, 7X, 

1 6HX "AXIS, 8 X , ft H V A R I A NO E , 3X, 6HY AXIS , 8X , 8H V ART ANC £ , 8X , 

L 6HZ AXIS, ftX, 8H VARIANCE,/) 

C CALCULATE MOMENTS 

RP - R ( 1 ) * ’ ( 1 ) *R ( 1 ) * . 0 0 0 0 1 
ALPH A (1 ) ~ A L PH A ( 1) *RR 
BETA (1)-- BETA (i)*RR 
GAMMA ( i ) - -GAMMA (1 ) * R R 
I R = NE QX -N T E p MS 
DA (1) = DA ( 1 ) * RR* T ( I R) 

D P < 1 ) = 0 3 (1_) *RR* U IR1 _ 

I P- ML 07 -N T ERM2 

QC (1 ) =QC( 1 ) »RR»T (IR) 

WRITE (61, 20 50) ALPHA { 1) ,0 4(1) , BETA (1 ) ,OB (1) , GAMM A (1 ) , DC ( 1 > 

2060 FORMAT ( 4 ( / )_, 10 X, 2 bHC ALCULA TEO DIPOLE MOMENTS ,20X, 29H90 PERCENT 
1 COMF1 D ANCE L I M I T S ,7,2 OX, 9 H M NO RT H = , F2 0 . 10 , 20 X , F20 . 1 0 , / , 2 0 X , 

19HNEAST = ,F2G .10 ,?GX ,F20 .1 0 ,/ , 20 X , 9HM0 0 WN = , F 2 0 • 1 0 , 2 0 X , F2 0 . 

110,////) 

C CHECK FOR ANOTHER ANALYSTS __ 

~ N1FRMS = N TERMS +1 

NTFRM? = NTFRM2 + 1 
IF ( N E Q 7 - NTFRM 2) 3008,300 3, 3009 
30 0 8 NT E PM 2 = 1 

3009 IF (ME QX - NT ER MS ) 3010,301 0,411 

30 IP NTEPMS = 1 

GO TO 423 ~~ ~ ~~ ~ ' ' ’ ’ 
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43P CONTINUE 
STOP 

END 


5 J iRPnil T TNF HAT RI X (NCON , N A X I.5., NRj, N E-Q.V -IP-t > 

C TOFCRN WEIGHTED MATRIX EQUATIONS 

DIMENSION R(1),A(1), 0(1) __ 

DIMENSION C(25,25) , M(25,25) 

DIMEN SI O N Y ( 25 ) , V V ( 2 5) , X< 2 5 ) * R AD ( 4 ) 

DIMENSION Q (25,25) 

_ REAL M - - - 

DO 1 I = 1,25 

X(I) = Y(I) ~ VY(I) = 0, 

DO 1 J = 1,25 

1 C (I, J) = H(I , J) = 0. 

DO ? I = 1,NR 

2 RAD ( I ) = F(1)/R(I> - 

IF(NCON) 6,6 ,3 


3 T F (N AX IS ) 5,5,4 


4 

N = NEQ/? 
GO TO 9 




5 

N = NFO 
GO TO 9 




6 

7 

IF (NAXIS) 8,6,7 
N = ( NEQ- 1) 72 




9 

GO TO 9 
N = NEQ-t 




0 

DO 17 I = 1 , N 
DO 17 J = 1,IP 




RP = RAD { I ) ** ( 2* J +1 ) 
F = i. 


OF = 1. 

TJ - 1 . 

DO 10 K - 1,J 

F = F*F LO AT (K) 

D F = DF*FLOAT(2*K-l) 

10 TJ = - 2 » * T J 

FF = (-2* *FLO AT (J) ) / ( T J* F ) 

RR = RR*FF*OF 

IF (MAXIS) 12,12,11 

11 CC I , J) = 2. AFLOAT (J ) *RR 
I N = I + N 

C ( I N , J ) = RR 

GO TO 13 

12 C ( I , J ) - RR_ _ ...... 

13 CONTINUE 

I F ( NC ON) 14,14,17 

1 4 RA = RAD ( NR ) ** (2* J+l ) 

C COMPUTE C (NFQ,J) FOR U CONFIGURATION 
IFTNAXIS* 16,16,16 

15 C ( N E Q , J ) = FLOAT (2* -J -l) »FL OAT <J_>»RA 

” GO TO 17 

16 C (NEQ , J) = 2 « ♦ FLOAT C J^RA 

17 CONTINUE 
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C HAVE (C(I,J) FO R ALL CQNFIGS AND AXES 

~C CONSTRUCT LLSQ AND DATA MATRICES ' 

no is j = i,ip 

0 0 18 I = 1,NEQ " 

18 X < J ) - X < J ) + ( A (I ) *C_( I , J) )/(0(I) *0_( I ) > 

DO 19 K = i,IP ~ " ' " ' 

DO 19 J - 1,IP 

DO 19 I -i ( NEQ. — 

19 M(J,K> = M(J,K) 4 <C<I, J) *C( I*K> > /<D(I> *D< I) > 

C HAVE WEIGHT ED MATRICES ““ - 

DO HO I ' 1,IP 

Y(I) = X ( I) ' 

30 CQNTTN UF 

CALL MATTNV <M, IP, Y , i, DETER) 

JC Q IS NOW THE INVERSE OF MO) ( Y = MULTIPOLE MOMENTS 

C COMPUTE DEVIATIONS AND RESIDUALS 

RES = 0. 

SS2 = 0, 

SR - Qj. 

SQ = 0. 

DO 21 I = 1,NF Q 

FL = 0 • 

DO 22 J - It IP 

FL = FL * C(I,J)*Y(J) 

22 CON TINUE 

S = (A <I> -FL) /0( I) 

R FS = RES ± S*S 

SS = S/AC I) 

SS = SORT (SS»SS) 

SR = SR + S 
SO = SO f SS 
SS2 = SS? + SS*SS 

?1 CONTINUE 

IDF = NEC — I p 
S 2 ~ RES/ FLOAT (IDF) 

4L4 FORMAT ( 1 0 X , 6 ( E 1 5 . 7 , 2 X ) ) 

498 FORMAT ( 5 ) X , 12H0 A f A MATRI X , / ) 

491 FORMAT < 2X , / ,50 X, 20HMEA SUREMEN T MATRIX ,/) 

492 FORMAT ( ? X,/,50X,8HINV£RS£ , /) 

493 FORMAT ( 1 OX , 11H PARAMETERS' ,/> “ 

434 FORMAT ( IPX, 21HQEGRE ES OF FREEDOM = , 19,/, 

1 IPX, 16 4 SIGMA SQUARED = , E20.10 , /, 10X, 8H0ETM = , 

1 £ 20.10 , /) 

412 FORMAT < 1 1 X , / , 5 OX , 2 0 H NE AR FIELD MATRIX ~ ,//) 

C VARIANCE OF Y(I) = VY(I ) 

DO 23 I - 1 , IP 
QU = M (I, I) *S2 
I F ( OU • L T • 0 • ) 3 01,902 

9 11 WRITE (61, 90 3) 1,1, 1, M (I , I ) , RESUME 0 , 1 P , S 2 

503 FORMAT (IPX, 4 HD EV ( ,I2,1X,32H) IS NEGATIVE, PARAMETERS ARE..",/, 

1 1 0 X , 8 H I N V E RS E ( ,I ?> 1 H , ,12, 4H) = , F2 0 . 1 0 , / , 1 OX , 5HRES - , FRO. 14 

1 , / , 9 H N t 0 = , 1 2 , 2 X , 9 H i P = , 1 2 , ? X , 5 H S 2 = , F 2U . 14) 

QU = -OU 

9)2 VY(T) = SQRT(QU) 

23 CONTINUE 

RETURN "" ™ 

END 
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SUBROUTINE M AT IN V ( A , N, 8 , H , DETERM) 

DIMENSION I PIVOT (25) , A ( 25 , 25 ) , B ( 25 , 1 ) , I NOE X ( 25 , 2) 

C E 

C INITIALIZATION F 

C F 

□ETERM= 1, 0 

DO 20 J=1,N F 

2 0 IPIVOT ( J) =0 F 

DO 55 0 1=1, N F 

C F 

C SEARCH FOR PIVOT ELE MENT __ _ . F 

C F 

AN AX-0 > 0 F 

DO i Q 5 J-1,N 

IF ( IPIVOT (J)) 10 5,60 

60 DO 100 K“ 1 , N F 

IF (IPIVOT (K)-l) 60, 100, 740 

80 IF(A8SF (A (J,K> ) . GT . A BS F ( AM AX ) ) 85,100 

85 IPOW=J F 

ICOLUM=< 

A MA X = A ( J, K) 

100 CONTINUE F 

105 CONTINUE 

IPIVOT (10 OLUM) = IPIVOT (ICOLUM) +1 F 

C __ F 

C INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL F 

C E 

IF (I ROW, EQ. ICOLUM) 260,140 

V+n OFT FRM = -D ET CRM J 

DO 200 L = I, N 

SWAP = A (IROW,L) 

A ( IPOW, L) = A ( ICOLUM, L) 

200 A (ICOLUM, L) -SWAP F 

IF <M) 210 ,260 

21 0 DO 250 L-l, M 1 

SNA P = 8 ( IR OW , L) F 

B (IRON, L) =3 (ICOLUM, L ) _ F 

250 3 ( ICOLUM, L> =S~WAP F 

260 INDEX (1,1 )=IROW F 

INDEX (I ,2 )= ICOLUM F 

PIVOT = A< ICOLUM, ICOLUM ) 

DETFRM=OL TERM*PI VOT 
PIVTN- 1,0/ PIVOT 

C DI VIDE PI VOT ROW BY PI VOT ELEMENT “ " " “ F 

c f 

A (ICOLUM, ICOLUM) =1.0 

DO 350 L = 1 * N F 

36 0 A ( IC O L U M ,L)=A (ICOLUM, L)^PIVIN “ " “ 

IF (M) 360,380 

36 0 DO 370 U 1 » M ” ~ ' “ F 

3/0 B ( ICOLUM, L) =B( ICOLUM, L> * PI VI N 

C F 

C REDUCE NON- PIVOT ROWS _ F 

3 8 0 DO 550 LI _=lj_N F 

IF (LT'VEO .ICOLUM) 550,400 
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400 T = A (Lit 10 OLUM) 

A <L1 ,ICOLUM> =0 .0 

DO 450 L- 1 , N 

450 A < L 1 , L > =A < LI , L ) - A < ICTON 7LT*T“ 
IF (M) 46- ,550 
480 DO 500 L=i,H 


5>)Q 9 (LI, L ) =B_(JLlj ! .L)rB (ICOLUM, L ) *T 

580 CONTINUE 
C 


c 

_c_ 

INTERCHANGE COLUMNS 

F 

F 


DO 710 1=1, N 
L=N+1-I 

F 


IF (INDEX (L,l) .EQ.INDEX(L,2) > 710,630 
630 Jf?OW = INDF X (L,l) 

F 


JCOLUM= INDEX (L,2) 
DO 705 K=i,N 

E 


SWAP=A (K, JROW) 

A (K , JROW) =A (K, JCOLUM) 

F 


A (K, JCOLUM) =SWAP 
705 CONTINUE 

E 


710 CONTINUE 
740 RETURN 

F 

F 


END “ 



SUBJRfiUTlNE FOURIER ( NFIT • DATA. CO~F) 


c 

c 

CALCULATE FOURIER COEF UP TO AMPLITUDES OF 

SIN(NFIT*THETA) 

c 

c 

AND COS( NFIT*T HE TA ) 
COEF ( 1) =1 C TERM 


c 

c 

C0EFC2K) IS COEF OF COS (.<* THETA ) 
COEF ( 2K+1 ) IS COEF OF SIN ( K*THETA) 


c 

DIMENSION DATA <l),COEF(l) 


A - U . 

DO 5 1=1,36 

5 A=A+DATA(I) 

COEF (1) =A/36. 


OO 100 N= 1 , NFIT 
THETA=-.l 745327 


B1-T1=62=T2=0. 

DO 50 J = 1 ,36 


THETA=THh TA+. 1745327 
C = COS (N*7 HETA) 



S = S I N (N*T BETA) 
B?=B2+S*S 


o 1= H1 + C*C 

T1 =D AT A ( 1 )*C+T1 


->0 T2 = DATA (J) *S+T 2 

ISTOR=2*N 


COLF(ISTDR) =T1/B1 
too COEF ( I S TO R+ 1 ) = T 2 / □ ? 



END 
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lx. lx Ix lx 


SUBROUTINE LIST (AMP«NTERM» IFORMAT) 

OIMENSION AMP(i) 

GO TO ( 1 +2 *3) IFORMAT 

1 ID=1HX 

GO TO 4 

2 I Q= 1H Y 

GO TO 4 

3 I D= 1HZ 

4 J-2 

WRITE (61, 300) 10 , AMP ( 1 ) 

300 FORMAT (31X,2HDC, At, 3H = , E 4 ) 

00 400 1=1* NTERM 

WRITE (61 . 10 0) I, ID, AMP( J) , I , ID, AMP< J+i) 

100 FOR HAT <3nx*I2*lHC*Al*2H= , El 1 . 4 , 5X , I 2 , 1HS , A 1 , 2H= 

400 J=J + 2 

WRITE (61* 500) 

500 FORMAT (/) 

RFTURN 

END 


Ell. 4) 


SUBROUTINE MAXMIN ( AM AX , AMI N . N . PT S) 

OIMENSION PTS(l) 

AMAX=AMIN = PTS (1) 

DO 100 1=2, N 

AMAX=AMAX 1 ( AMAX , PTS (I) ) 

100 AMIN=AMIN1 (AMIN, PTS ( I) ) 

RETURN 

END 


NASA-Langley, 1972 
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